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SUMMARY 


The  purpose  of  the  present  report  is  to  investigate  the  hodograph 
method  as  it  is  applied  in  general  to  the  problem  of  compressible  flow. 
First,  the  hodograph  equations  are  given  in  various  canonical  forms 
which  are  convenient  for  obtaining  solutions  in  the  different  flow 
regimes. 

Since  the  coefficients  of  the  canonical  differential  equations  are 
implicit  functions,  exact  solutions  are  difficult  to  find.  Consequently 
different  approximations  are  chosen  so  that  some  simpler  differential 
equations  capable  of  solution  can  be  obtained.  For  most  of  the  cases, 
fundamental  or  singular  solutions  are  given  or  indicated. 

The  detailed  development  is  concentrated  on  Chaplygin' s  second 
equation.  The  first-order  approximation  is  well-known  as  the  Tricomi 
equation.  The  second-  and  third-order  approximations  have  a  rather  new 
approach.  Both  approximations  follow  the  exact  gas  law  closely  in  the 
neighborhood  of  the  sonic  velocity.  The  solutions  are  found  to  be 
Whittaker  functions  and  the  associated  confluent  hypergeometrical 
functions.  Both  approximations  can  be  applied  to  the  incompressible 
flow  so  that  Chaplygin's  procedure  of  borrowing  the  boundary  conditions 
can  be  used  if  necessary.  For  the  third-order  approximation,  the  corre¬ 
sponding  hypothetical  gas  law  is  derived  and  is  found  to  differ  very 
little  from  the  exact  gas  law.  The  transformation  relation  between  the 
hodograph  plane  and  the  physical  plane  is  also  given  for  the  various 
solutions  considered. 

To  make  a  comparison  of  the  present  approximate  solution  with  the 
exact  Chaplygin  solutions,  the  flow  through  an  aperture,  studied  by 
Chaplygin  and  Lighthill,  is  reexamined.  There  is  some  difference  in 
the  problem  itself,  as  well  as  in  the  method  of  Chaplygin  and  Lighthill. 
First,  the  vessel  with  straight  walls  inclined  at  an  arbitrary  angle  is 
considered  rather  than  that  with  the  wall  at  right  angles.  Second,  no 
association  of  the  boundary  conditions  with  those  for  incompressible 
flow  is  made.  The  problem  is  treated  directly  as  a  boundary-value 
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problem.  The  result  calculated  with  the  Whittaker  function  checks 
well  with  that  obtained  by  Chaplygin  and  Lighthill. 


1  -  INTRODUCTION 


For  the  isentropic  steady  two-dimensional  flow  of  nonviscous  com¬ 
pressible  fluids,  the  pair  of  differential  equations  in  terms  of  the 
stream  function  and  the  potential  are  nonlinear ,  and  cannot 
be  solved  analytically  except  for  a  few  special  cases.  Chaplygin 
(reference  l)1  introduced  the  hodograph  method  in  order  to  transform 
these  nonlinear  equations  to  linear  ones,  so  that  the  available 
classical  mathematical  analysis  and  the  principle  of  superposition  can 
be  applied.  While  there  is  a  gain  in  the  linearity  of  the  equations, 
some  new  difficulties  arise  in  the  hodograph  method.  First  of  all, 
it  is  difficult,  in  general,  to  transfer  the  physical  boundary  condi¬ 
tions  to  the  hodograph  plane.  Second,  the  flow  in  the  hodograph  plane 
usually  consists  of  multiple— sheeted  Riemann  surfaces  with  a  number 
of  singularities  where  the  analytic  continuation  of  the  series  solu¬ 
tions  becomes  very  complicated.  Third,  only  one  kind  of  a  particular 
solution  of  the  stream  function  -  hypergeometrical  functions  and 
trigonometrical  functions  -  has  been  obtained  so  far  for  the  equations 
in  the  hodograph  plane.  By  superposition,  a  series  solution  is 
achieved,  but  is  difficult  to  apply  even  if  the  boundary  conditions 
are  known.  Fourth,  the  transformation  between  the  physical  and  hodo¬ 
graph  planes  becomes  singular  when  the  Jacobian  determinate  becomes 
zero  or  infinite.  This  is  likely  to  happen  when  the  supersonic  flow 
is  imbedded  in  a  subsonic  flow  region.  The  existence  of  the  so-called 
limiting  line  (reference  7)  causes  a  breakdown  of  the  entire  flow. 

In  order  to  obtain  some  solutions  in  the  case  of  subsonic  flow, 
Chaplygin  discovered  an  ingenious  method  of  obtaining  usable  solutions 
for  compressible  flows  by  comparison  with  the  corresponding  series 
solutions  of  the  incompressible-flow  patterns  as  the  limiting  case  in 
the  hodograph  plane.  In  the  case  of  no  circulation  around  the  closed 
body,  the  derived  series  solutions  involve  a  free  constant  to  be  chosen 
at  pleasure.  But  fortunately  each  solution  chosen  corresponds  to  a 
reasonable  flow  pattern  in  the  physical  plane.  As  the  free- stream 
Mach  number  increases,  the  body  shape  deviates  from  the  image  body  in 
the  incompressible  flow.  Thus,  in  general,  the  body  shape  cannot  be 
preassigned  and  the  compressible  flow  about  it  determined  at  a  certain 
Mach  number.  This  method  has  been  carried  further  by  Tsien  and  Kuo 
(references  2  and  3),  Lighthill  (reference  4),  and  Cherry  (reference  5) 
on  transonic  flow  up  to  the  occurrence  of  the  limiting  line.  Some  of 
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their  solutions  apply  to  simple  closed  "bodies  with  circulation.  There 
are  difficulties  in  the  analytic  continuation  for  an  arbitrary  flow 
pattern. 

Although  Chaplygin  and  others  were  pioneers  in  developing  and 
adapting  the  hodograph  method  to  gas  dynamics  in  the  early  part  of  this 
century,  results  of  the  later  Russian  investigators  were  not  available 
in  the  English  language  up  to  19*40 .  Since  the  second  World  War,  the 
work  of  Russians  has  again  assumed  prominence  in  the  field,  namely, 
Frankl  (reference  ll),  Falkovich  (reference  3*0,  and  Christianovitch 
(reference  *4-2) .  Along  the  line  of  formulation  by  Tricomi,  Frankl' s 
proof  of  the  existence  and  the  uniqueness  of  the  solution  of  Chaplygin's 
equation  in  the  transonic  regime  is  interesting.  Frankl  and  Falkovich 
obtained  the  solution  for  the  channel  flow  and  showed  that  the  stream 
function  is  triple-valued  in  the  hodograph  plane  of  the  axially  sym¬ 
metrical  plane  at  sonic  velocity.  These  results  confirm  what  has  been 
achieved  by  Lighthill  (reference  4).  Recently,  Tomotika  and  Tamada 
(reference  33)  have  formulated  some  approximate  nonlinear  hodograph 
equations.  A  number  of  interesting  particular  solutions  for  the 
channel  flow  have  been  obtained.  Of  course,  the  solutions  are  not 
superimposable.  Ehlers  (reference  10)  and  Carrier,  along  the  line  of 
Christianovitch,  have  obtained  the  fundamental  solutions  of  the  Tricomi 
equation  and  the  corresponding  channel  flow  to  the  second  order  of 
approximation. 

With  the  notion  of  the  correspondence  of  the  incompressible  flow 
to  the  compressible,  Bergman  (reference  6)  has  developed  an  integration 
method  for  calculating  the  subsonic  flow.  For  the  supersonic  case  an 
extension  of  the  Riemann  method  is  also  made  by  an  iteration  process 
(reference  4l).  Bers  and  Gelbart  (reference  55)  have  similarly  developed 

a  line-integral  operator  to  construct  the  so-called  r  -monogenic  com¬ 
plex  function  which  satisfies  the  hodograph  equations  or  the  general 
Cauchy- Riemann  equations  in  equivalence.  All  are  important  contribu¬ 
tions  to  the  solution  of  the  differential  equation  but  offer  rather 
difficult  ways  to  obtain  useful  solutions  for  the  flow  of  compressible 
fluids . 

In  the  last  few  years,  Guderley  (reference  30)  and  Yoshihara 
(reference  36)  have  given  a  number  of  papers  on  transonic  flows, 
particularly  with  the  application  of  the  transonic  similarity  law  which 
was  developed  independently  by  Von  Karman  (reference  29)  and  Guderley 
(reference  30) .  They  have  achieved  very  important  approximate  results 
in  transonic  flow. 

In  the  subsonic  case,  if  a  linear  approximation  is  used  to  replace 
the  isentropic  pres sure- density  relation,  Karman  (reference  7)  and 
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Tsien  (reference  l4) ,  following  another  approach  of  Chaplygin,  have ' 
obtained  many  useful  results,  particularly  the  well-known  K^rman-Tsien 
formula  (reference  7)  for  pressure  corrections.  Lin  (reference  15)  has 
discussed  the  conditions  for  obtaining  closed  contours  in  the  physical 
plane,  if  circulation  exists.  Clauser  (reference  2k)  recently  applied 
this  method  to  find  the  body  shape  with  preassigned  pressure  distribu¬ 
tion  and  this  is  definitely  very  useful  in  high-subsonic  laminar-flow 
problems.  Garrick  and  Kaplan  (reference  20 )  have  another  interesting 
approach  to  pressure-correction  formulas  using  Ringleb's  solution 
(reference  8)  of  a  simple  source  and  vortex. 

The  present  investigation,  as  the  initial  step  of  the  research 
program  of  this  challenging  problem,  is  mainly  interested  in  the 
following  three  aspects  of  the  problem: 

(a)  Besides  Chaplygin's  differential  equations  (reference  l), 
can  other  useful  forms  of  the  differential  equations  be  found  system¬ 
atically,  particularly  the  canonical  forms  in  the  different  flow  regimes? 
Of  particular  interest  are  the  fundamental  or  singular  solutions  which 
represent  the  types  of  singularities  encountered  in  the  hodograph  plane. 

Is  there  any  method  of  constructing  such  solutions  as  shown  by  Picard, 
Hadamard  (reference  52),  Hilbert,  Riemann,  and  Tricomi  (reference  53)? 

(b)  If  such  fundamental  solutions  are  too  difficult  to  construct 
or  too  complicated  to  apply  to  the  flow  problem,  what  other  reasonable 
approximations,  besides  the  Karman-Tsien  (reference  7)  approach,  can 
be  made  so  that  some  useful  results  can  be  derived  in  the  different 
flow  regimes? 

(c)  It  is  well-known  that  the  boundary  conditions  of  the  stream 
function  (not  the  potential  0)  are  well-defined  in  the  hodograph 
for  polygonal  profiles  with  or  without  free  streamlines  up  to  transonic 
flow.  As  a  reasonable  preliminary  approach,  can  the  compressible  flow 
around  such  a  given  body  be  found? 

In  brief,  this  paper  contains  a  systematic  list  of  useful  forms  of 
the  differential  equations.  The  canonical  forms  of  the  exact  differential 
equations  are  given  for  subsonic,  supersonic,  and  transonic  regimes. 
Unfortunately,  one  of  the  coefficients  of  the  equation  is  an  implicit 
function  of  one  independent  variable.  This  makes  it  impossible  to  con¬ 
struct  fundamental  solutions  which  would  be  of  value  for  practical 
application. 

Approximations  to  the  implicit  function  are  chosen  in  such  a  wny 
that  the  solutions  can  be  found  from  classical  mathematical  analysis. 

In  applying  the  approximations,  three  objectives  are  kept  in  mind: 

First,  the  differential  equation  must  reduce  to  the  Laplace  equation 
as  M  — >  0  or  M  - >  M,,,  so  that  Chaplygin's  procedure  of  utilizing 
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the  incompressible  flow  can  be  followed.  Second,  in  the  transonic 
range,  the  transonic  similarity  law  of  Karman  (reference  29)  can  be 
applied  and  the  simplified  boundary  conditions  can  be  used.  Third 
the  singular  solutions  and  the  fundamental  solutions  in  the  subsonic 
and  transonic  regime  can  be  found  so  that  the  flow  can  be  determined 
directly  if  the  boundary  conditions  in  the  hodograph  plane  are  assigned. 
Or,  alternately,  the  Riemann  function  in  the  supersonic  regime  can  he 
found,  so  that  an  integral  solution  can  be  obtained  when  the  initial 
value  or  the  Cauchy  data  along  a  noncharacteristic  line  are  specified. 
Last,  the  characteristics  method  can  he  used. 

In  the  present  paper,  a  number  of  approximations  for  the  canonical 
forms  in  the  subsonic  flow^are  given.  The  zero-order  approximation  is 
actually  the  same  as  Von  Karman's  approximation  given  in  reference  7, 
equation  (63),  or  reference  13,  pages  186  to  188.  Both  differential 
equations  of  the  first-order  approximation  can  be  reduced  to  the 

Laplace  equation  in  polar  coordinates.  The  singular  and  fundamental 
solutions  have  heen  given. 

The  second-order  approximations  to  the  differential  equations  can 
he  reduced  to  Stratton*  s  equation  (reference  70)  hy  means  of  the 
separation  of  variables.  The  particular  solutions  are  indicated. 

There  are  some  better  approximations  which  should  hold  for  any 
subsonic  Mach  number.  They  are  shown  in  section  3,  "Solutions  to 
Canonical  Forms  of  Approximate  Differential  Equations  in  Subsonic, 
Supersonic,  and  Transonic  Regimes."  Similar  approximations  are  obtained 
m  the  supersonic  region.  Of  course,  for  each  approximation,  the 
corresponding  pressure-density  relation  must  be  determined. 

If  the  general  class  of  singular  solutions  in  the  hodograph  plane 
can  be  obtained,  then  it  should  be  possible  to  solve  a  large  number  of 
problems  by  a  suitable  placement  of  the  singularities  in  the  hodograph 
plane  to  satisfy  the  desired  boundary  condition  in  much  the  same  way 

that  sources,  vortices,  doublets,  and  so  forth  are  used  in  incompres- 
sible  flow,  * 


Some  investigations  on  possible  approximations  to  Chaplygin's 
second  equation  are  made  in  section  4,  "Different  Approximations  to 
Chaplygin  s  Differential  Equation  and  Their  Solutions."  This  second 
equation  is  more  convenient  to  use  in  flow  than  the  first  equation 

rt;°rder  ?Pproximation  to  this  differential  equation  has  been 
shown  by  Frankl  (reference  11 )  to  be  Tricomi's  equation.  Both  the 
second-  and  third-order  approximations  are  shown  to  be  associated  with 
Whittaker  s  equation  (reference  6l).  The  third  approximation  is  suf- 
ficiently  good  that  it  should  give  fairly  accurate  results  in  all 
transonic-flow  problems.  The  corresponding  hypothetical  gas  law  is 
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shown  in  section  5,  "Hypothetical  Gas  Law  Corresponding  to  Approximations 
of  Chaplygin's  Second  Equation,"  and  deviates  so  little  from  the  exact  gas 
law  in  transonic  range  that  the  exact  gas  law  may  justifiably  replace 
it.  The  transformation  relation  between  the  physical  and  hodograph 
planes  is  also  shown  for  the  third  approximation. 

The  polygonal  body,  either  closed  or  open,  is  known  to  give  simple 
boundary  conditions  in  the  hodograph  for  the  stream  function  Jr  but 
not  for  the  potential  0.  The  particular  solutions  obtained  by  the 
principle  of  separation  of  variables  and  the  series  solutions  obtained 
by  superposition  cannot  be  applied  directly  without  solving  systems  of 
an  infinite  number  of  simultaneous  equations  to  determine  the  infinite 
number  of  coefficients  in  the  series.  In  this  class  of  bodies  with 
one  or  more  convex  corners,  the  solution  of  the  incompressible  flow 
cannot  be  used  because  the  velocity  at  such  a  corner  is  infinite 
according  to  the  theory  of  incompressible  flow.  Physically,  the  flow 
passing  such  a  body  is  always  transonic  in  character,  no  matter  how 
low  the  free- stream  velocity  is.  As  shown  by  Guderley  (reference  31) 
and  Busemann  (reference  38)  some  shock  always  occurs  at  such  corners. 

But  with  the  simplified  assumption  of  the  transonic  similarity  law, 
the  flow  about  this  group  of  bodies  should  be  obtainable.  In 
section  3,  "Solutions  to  Canonical  Forms  of  Approximate  Differential 
Equations  in  Subsonic,  Supersonic,  and  Transonic  Regimes,  a  more  gen¬ 
eral  simplified  equation  satisfying  the  transonic  similarity  law  is 
given  and  the  solutions  are  shown. 

For  an  open  body  built  with  straight-line  elements  but  with  no. 
convex  corners  the  problem  can  be  attacked  with  the  Chaplygin  technique 
or  solved  directly  as  a  boundary- value  problem.  In  this  case,  no  free 
constants  can  be  chosen,  and  the  solution  is  uniquely  determined. 

Owing  to  the  complicated  nature  of  the  asymptotic  behavior  of  the 
Whittaker  function  with  a  very  large  parameter,  and  simultaneously  with 
a  very  large  value  of  the  independent  variable,  the  future  work  will 
devote  a  considerable  amount  of  time  to  finding  the  asymptotic  solu¬ 
tions  corresponding  physically  to  different  flow  regimes.  This  is  an 
important  step  in  solving  flow  problems  with  the  Whittaker  function. 

The  Whittaker  function  converges  very  slowly  for  a  large  parameter,  and 
is  quite  similar  to  the  Chaplygin  function  in  this  respect. 

In  order  to  demonstrate  the  method,  the  flow  through  an  aperture 
is  reexamined  in  order  to  check  with  the  results  of  Chaplygin  and 
Lighthill.  It  is  found  that  the  third-order  approximation  agrees  very 
well  with  their  work.  This  is  given  in  the  last  section  of  the  paper. 

This  investigation,  conducted  at  The  Johns  Hopkins  University,  was 
sponsored  by  and  conducted  with  the  financial  assistance  of  the  National 
Advisory  Committee  for  Aeronautics. 
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2  -  CANONICAL  AND  OTHER  FORMS  OF  DIFFERENTIAL  EQUATIONS  IN 
SUBSONIC,  TRANSONIC,  AND  SUPERSONIC  REGIMES 
2.1  -  General  Transformation  of  Differential  Equations  in 

Hodograph  Plane 


The  pair  of  Chaplygin's  first-order  simultaneous  partial  differ¬ 
ential  equations  for  the  stream  function  i  and  the  potential  0  in 
the  hodograph  plane  (equations  (11.10)  and  (11. 11)  in  reference  13)  is 


=  _  Po  1  -  M2  M 
dq  p  q  .  be 


=  Po 

80  p  ^  Sq 


(lb) 


where  pQ/ p  and 

’  P 

variable  q. 


1  -  M2 


are  given  functions  of  the  independent 


The  above  system  can  be  transformed  to  many  other  forms  which  are 
more  adaptable  for  analysis  if  a  new  independent  variable  Q  is 
introduced  to  replace  q,  thus: 


Q  =  Q(q) 


list  of  symbols  is  given  in  appendix  A. 


(2) 
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In  terms  of  Q,  equation  (2)  yields 


80  dQ  _  Po  1  -  M2  8t 
8Q  dq  P  q  80 


(3a) 


_  £0  q  dQ  8t  (3b) 

8©  p  dq  8Q 

Eliminating  0  and  then  t  from  the  above,  a  pair  of  second-order 
partial  differential  equations  is  obtained  as  follows: 


82j  1  -  l^/dQ\~2  82j  8t  _d_ 

8q2  +  q2  W  8e2  e 


0 


(4a) 


1  -  M2/dQ\-2  1.  <  ^  fo  1  -  ^fjgV1  .  =  0  (4b) 

q2  \dq)  a02  8Q  dQ  e  Lp  <1  W /.  J 


Of  course,  if  Q  is  chosen  equal  to  q,  the  following  pair  of  equations 
is  obtained: 

8^  +  1  -  M2  82t  +  1  +  m2  _  0  ( ) 

8q2  q2  802  <1  ^ 


820  ,  1  -  M2  820  d  /pQ  1  -  m2\80  _  0  (ltd) 

8q2  q2  80 2  d(l\P  <1  M 


which  have  been  used  by  Ringleb  (reference  8).  Now  make  some  proper 
choice  of  the  function  Q(q),  so  that  either  one  or  both  of  the  above 
equations  can  be  transformed  to  some  simpler  form. 
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2.2  -  Chaplygin's  Differential  Equations  -  First  Form 

In  Chaplygin's  differential  equations  of  the  first  form,  he 
introduced  the  algebraic  relation  between  Q  and  q  as 


Q(q)  =  t  = 


7-1  q‘ 


(ao*) 


where  aQ*  is  the  stagnation  sound  velocity.  It  is  obvious  that  the 

range  of  t  ( t  has  an  important  physical  meaning.  It  is  the  ratio 
of  the  kinetic  energy  of  the  gas  to  the  total  energy  or  enthalpy)  is 

from  0  to  1.  Figure  1  shows  the  relation  of  t,  — —  ,  and  Mach  num- 

a  * 
o 

ber  M,  Equation  (4a)  transforms  to  his  well-knovn  equation 

a  f  2t  afl  ,  1  -  ^it  a2t  „  ,, 

3  8t  p+l  ^  2  1 

L(1  -  t)  J  2t(1  -  t)  00 


where  6  =  - . 

7-1 

i  m2  _  1  ”  ^1T 


7+1  r 

=  - ,  and  the  relations  — 

7-1  P, 


(l  -  t)P  and 


p  JL  —  |JL-l  \ 

1  -  M  =  - —  are  introduced.  By  means  of  the  principle  of  the 

1  -  T 

separation  of  variables,  the  particular  solutions  of  this  equation  are 
combinations  of  hypergeometric  and  trigonometric  functions.  The  dif¬ 
ferential  equation  in  terms  of  the  potential  0  transforms  to  a  much 
more  complicated  form 


a_  2t( i  -  t)p+1  b[\  +  (i  -  T)p  a20 
1  -  p-jT  8t  2t  £q2 


and  its  solution  is  difficult  to  obtain  directly,  but  it  can  be  obtained 
by  integrating  the  differentials  of  i|r.  Actually,  there  are  four  sin¬ 
gularities  of  this  equation  -  three  regular  ones  at  t  =  0, 

t  =  —  =  and  r  =  1  and  one  irregular  singularity  at  t  = 

Pj  23  +  1 

Therefore  its  solution  is  not  so  simple  as  the  other  equation. 


00, 
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2.3  -  Chaplygin’s  Differential  Equations  -  Second  Form 
To  eliminate  the  last  term  in  equation  (4a),  choose 


Q(q)  =  o 


(7) 


where  a  is  defined  by  (cCT  is  a  constant) 


(8) 


or 


P  P  dq 


Substituting  into  equation  (4a),  the  last  term  drops  out  and 


(9) 


_K_ 

where 


( 10a) 


(11) 


In  figure  2,  K  is  given  as  a  function  of  t .  Of  course,  K  should 
be  expressed  in  terms  of  d.  Unfortunately,  it  is  an  implicit  func¬ 
tion  of  a  and  the  differential  equation  is  impossible  to  solve 
exactly.  Figure  3  shows  the  behavior  of  K  as  a  function  of  a  for 
the  case  c^  =  -1  and  7  =  1.4  where  the  upper  limit  of  the  integra¬ 
tion  is  q  =  a*,  the  sound  velocity.  It  is  of  interest  to  observe  the 
two  asymptotes  of  the  K  -  o  curve:  (a)  K  — >■  1  as  M  *■  0  or 

cr  *co  and  (b)  K — >  -00  as  a  — >-0.2513-  The  function  a  can 

he  integrated  as  follows: 
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a 


Pa* 

I  P  <*q 


=  -0.2513  +  tanh_1( 1  -  t) 


1/2 


4- 


1  -  l  +  (1  -  t)2 

3  5 


(12) 

where  T  is  given  in  equation  (5)  and  7  =  1.4  is  explicitly  intro¬ 
duced,  In  figure  4,  a  is  given  as  a  function  of  t  .  If  equation  (9) 
is  substituted  into  equation  (4b), 


which  is 


rather  more  complicated. 


(10b) 


2.4  -  Another  Form  of  Chaplygin’s  Second  Differential  Equation 
To  eliminate  the  'last  term  in  equation  (4b),  introduce 


Q(q)  =  v 


(13) 


where  v  is  defined  by  (cv  is  a  constant) 


or 


Substituting  in  equation  (4b),  the  last  term  drops  out  and 
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where  K  is  the  same  as  defined  in  equation  (ll)  and  cv  =  -1  is 
chosen.  Of  course  K  should  he  expressed  in  terms  of  V,  if  it  were 
not  an  implicit  function  of  v.  No  detailed  investigation  will  he 
given  in  this  report. 

Similarly  it  can  he  shown  that 


K  &  +  tl  + 

8v2  8e2 


(16a) 


It  should  he  noted  that  equations  (10)  or  (l 6)  are  simpler  for  only  one 
of  the  pair;  at  the  sacrifice  of  the  other. 


2.5  -  Canonical  Forms  in  Subsonic  Flow 

It  is  well-known  that  the  fundamental  solution  may  be  found  when 
the  differential  equation  is  reduced  to  the  canonical  form.  Therefore, 
the  canonical  form  is  worth  while  to  obtain.  If  the  velocity  is  sub¬ 
sonic  everywhere  in  the  domain,  the  canonical  forms  of  the  differential 
equations  can  be  found  in  equation  (4)  in  the  subsonic  or  elliptic 
range.  Take 


Q(q)  =03 


(17) 


such  that,  for  M  <  1, 


or 


(i  -  ^)i/2 . 
dco  =  -  - - - -  dq 


( l8a) 


( Actually  do>  =  Here  the  minus  sign  is  chosen  so  that 

^  \  m 

a*  can  be  imposed  as  an  upper  limit. J  After  integration, 
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l£=  d(loge  K1/2) 

2  K  (to 

=  -(m  -  l)djT2 

(1  -  ^lT  )3^2( 1  '  T)1/2 

_  _J^1 _ .  (20) 

“1  -  1  (!  .  rf)3/2 

d(loge  K1/2) 

Figure  6  shows  — -  as  a  function  of  to. 

to 

Unfortunately  the  coefficient  of  the  first  derivative  is  an 
implicit  function  of  o>  and  also  is  singular  at  M  =  1.  Wiis  prevents 
the  use  of  Hilbert's  and  Hadamard's  approach  (reference  52)  to  formulate 
the  fundamental  solution.  Bergman  (reference  6),  Bers  (reference  IT); 
and  Gelbart  (reference  27)  in  this  country  and  Eichler  (reference  5») 
in  Germany  have  given  the  integral  solution  of  equation  (19)  but  the 
process  is  very  complicated.  The  present  paper  will  give  some  solu¬ 
tions  of  the  above  differential  equations  (19)  under  certain 
approximations. 

2.6  -  Canonical  Forms  in  Supersonic  Range 
If  the  velocity  q  in  the  flow  is  entirely  supersonic ^  take 


Q(q)  =  fi 


where  0,  is  defined  by 


M2  -  l/dn\2 


m - 


(22a) 


or,  if  the  positive  sign  is  chosen 


(m2  -  1)1/2 


d£2  = 
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In  the  hodograph  plane,  the  characteristics  are  the  same  as  those  of 
the  simple  wave  equation  and  are  invariant  to  the  behavior  of  K(ft). 

If  and  0  =  ( -K)^^0*  are  introduced,  the  following 

equations  are  obtained: 


5V 

afjj  +  Hie.)2  .  i  £  **  -  0 

(23c) 

5ft2 

520* 

ifh.)2  + 1  .  0 

(23d) 

5ft2 

Bq2  |_  l6\K  /  4  K_r 

a[ioge 

(-K)1^  -  l).^2 

( lilT  -  l)3/2U  "  T)V2 


(24) 

^1  "  1  (f^  -  l)3//2 


Figure  7  gives 


[loge  (-K)1/2] 


as  a  function  of  ft.  The  Riemann 


dft 


function  (reference  52)  always  exists  uniquely  and  the  solution  of  both 
equations  can  be  constructed  theoretically  for  a  region  bounded  by 
characteristics  if  the  Cauchy  data  are  given  along  a  line  that. is  not 
a  characteristic.  Owing  to  the  implicit  nature  of  the  coefficient 

i-ilog  (-K)1//21  as  a  function  of  ft,  the  Riemann  function  and  the 

dn  * —  ^  — 

solution  will  be  too  complicated  to  construct.  Besides  for  the  super¬ 
sonic  flow  the  Cauchy  data  in  the  transformed  hodograph  plane  are  not 
completely  known  in  general.  Therefore  the  Riemann  function  is  not 
of  much  use  in  the  hodograph  method. 


The  present  paper  will  also  show  some  solutions  to  the  differential 
equations  under  certain  approximations.  It  should  be  noted  that  ft  . 
and  cd  are  related  by  ft  =  Icd  if  it  is  desired  to  extend  the  defini¬ 
tion  of  cd  to  the  supersonic  side.  Such  an  extension  seems  very 
obscure  as  yet  in  its  meaning.  Fortunately  (dft)2  always  occurs  m 
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differentiation,  and  the  derivative  of  the  logarithmic  term  does  not 
depend  on  such  an  extension.  Therefore  equations  (23a)  and  (23b)  are 
single-valued. 

There  is  another  interesting  feature  of  the  canonical  form.  In 
the  hyperbolic  region,  the  characteristics  are  the  same  as  those  of  the 
simple  wave  equation  while  in  the  elliptic  region  the  imaginary  char¬ 
acteristics  are  the  same  as  those  of  the  Laplace  equation.  They  are 
invariant  to  physical  conditions.  The  physical  law  influences  only 
the  first-derivative  term. 


2.7  -  Canonical  Forms  in  Transonic  Domain 
Now  choose  Q(q)  =  6  where  e  is  defined  by 


(25) 


wt 


and  where  f(6)  is  an  arbitrary  function.  (The  upper  sign  is  to  be 
used  for  subsonic  flow  and  the  lower  sign  for  supersonic  flow. )  Then, 


[f(ej] 


1/2  r  [±  (i  -  m2  j] 


de  = 


dq 


For  M  <  1  choose 


f(e 


■©  ■ 


'\2  1  -  M2 


Equations  (4a)  and  (4b)  can  be  written  as 


(26) 


(25a) 


2L|+f(€.)8|+lLjL<log, 

6e'2  6e2  de'  de' 


1/2 


V  =  0 


(27a) 


(27b) 
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0  - f(e)  0  - 1  H-K)f(^01/2}  -  0  <®» 

which  are  valid  for  both  M  <  1  and  M  >  1.  It  should  he  noted  that 
e  >  0  for  M  >  1  and  e  <  0  for  M  <  1. 

The  canonical  form  of  the  mixed  differential  equation  is  explained 
as  follows.  If,  following  Tricomi  (reference  52)  f(e)  is  set  equal 

to  e  and  -K  =  ^  (m2  -  l)  is  introduced,  equation  (29)  reduces  to 
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a2* 

8e2 


1  at  _d_ 

2  de  de 


log, 


(30a) 


a H 

ae2 


a20 

ae2 


1  m  i.rog 

2  8e  de  L 


(30b) 


which  are  valid  for  the  mixed  region.  They  can  be  obtained  directly 
by  setting  ^  in  equation  (23)  if  M  >  1  and  ^  cuj 


in  equation  (19)  if  M  <  1.  These  may  be  called  the  canonical  forms 
of  the  differential  equations  in  the  transonic  region.  The  solution  of 
the  exact  differential  equation  is  obviously  difficult.  Some  singular 
solutions  of  the  above  equation  under  approximations  can  be  obtained 
as  will  be  shown  later. 


It  is  well-known  that  the  characteristics  are  fixed  for  all  linear 
hyperbolic  differential  equations.  But  there  is  an  interesting  feature 
in  the  canonical  form  such  as  equations  (30a)  and  (30b),  because  the 
characteristics  are  fixed  and  invariant  to  the  function  K(c)  or  to 
the  physical  problem  from  which  the  differential  equation  is  derived. 
The  characteristic  equations  are 


(e  -  Oo)2  -  j*3  -  0 

which  represent  two  families  of  characteristics  of  cubic  parabolas  with 
cusps  at  their  points  of  intersection  with  the  €  =  0  axis. 

Figure  8  shows  .6  as  a  function  of  go  and  Q.  Figure  9  shows 

€  |l°ge  as  a  function  of  e. 


3  -  SOLUTIONS  TO  CANONICAL  FORMS  OF  APPROXIMATE  DIFFERENTIAL 
EQUATIONS  IN  SUBSONIC,  SUPERSONIC,  AND  TRANSONIC  REGIMES 


In  the  last  section,  the  general  derivatives  were  shown  of  the 
canonical  forms  of  the  differential  equations  in  the  different  regimes. 
At  the  same  time,  it  is  found  that  the  coefficient  of  the  first- 
derivative  term  in  any  one  of  the  canonical  forms  is  always  an  implicit 
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function  of  one  independent  variable.  The  exact  solutions  are  very 
difficult  to  obtain.  It  seems  justifiable  to  seek  some  differential 
equations  with  an  approximate  coefficient  for  the  first-derivative 
term  so  that  solutions  are  relatively  easy  to  formulate.  Of  course , 
the  validity  of  the  solutions  thus  obtained  must  be  examined  and  the 
limits  of  the  range  of  application  must  be  defined.  In  this  section, 
the  solutions  of  the  approximate  differential  equations  will  be  given 
and  their  applications  will  be  discussed. 

3.1  -  Subsonic  Canonical  Forms  -  Equations  (19a) 

and  (19b) 

In  the  subsonic  case  three  approximations  to  the  canonical  form  of 
the  differential  equation  will  be  given.  Figure  10  shows  the  exact 

curve  of  =  ^.(l  -  ^  against  ox 

P 

Zero-order  approximation  of  Von  Karman.-  In  reference  Von  Karman 
suggests  taking 


K1/2  =  ^(l  -  M2)1/2  ~  ~(l  -  Ko2)1/2  .  (31) 

where  p  and  M^,  are  the  density  and  Mach  number  at  the  free  stream. 
In  figure  10,  the  horizontal  dotted  line  shows  the  approximation 

K1/2  .  It  is  apparent  that  the  approximation  becomes  better  as 

o>=cq* 

— >  0  or  M  — >  0.  With  this  approximation,  the  coefficient  of  the 
first  derivative  becomes  zero  in  equation  (19)  and  so 


d2t  +  dfi 

cto2  de2 


(32a) 


¥l+  =  o  (32b) 

5m2  Be2 


which  are  Laplace  equations.  These  equations  are  invariant  under 
translation  and  rotation.  Therefore,  <x>  can  be  set  equal  to  0  at 
sonic  velocity  a*  and  so  defined  as 
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CD 


P**  (l  -  M2)1^ 

Jq  q' 


(33) 


Of  course,  Chaplygin’s  procedure  of  utilizing  the  incompressible-flow 
pattern  to  find  the  corresponding  compressible  flow  can  be  followed. 
There  is  another  possible  approach.  If  the  types  of  singularities  and 
their  locations  in  the  flow  of  the  hodograph  plane  are  known  and  if, 
in  addition,  the  boundary  conditions  of  the  stream  function  are  known 
under  reasonable  assumption  or  under  assignment,  the  solution  of  the 
flow  should  be  obtained  uniquely  by  superposition  of  the  fundamental 
solutions  corresponding  to  the  types  and  location  of  the  singularities. 
After  the  solution  in  the  hodograph  plane  is  obtained,  the  corresponding 
flow  in  the  physical  plane  should  be  checked  to  determine  if  it  is  pos¬ 
sible  or  not.  To  serve  such  a  purpose,  the  fundamental  solutions  of  \|r 
and  0  due  to  a  source  at  o>0  and  Qq  are  given  as  an  example. 


0(co,6;  odo,0o)  =  -loge  jcn  -  cct,)2  +  (0  -  90)^2  (34a) 


/  \  CD  -  0) 

i(r(a),e;  <%,e0)=  -tan"1  - - —  (34b) 

e  -  0O 

Of  course,  with  the  simple  source,  there  can  be  built  up  the  potential 
and  the  stream  function  of  higher-order  singularities  such  as  doublets 
and  quadrupoles.  Actually,  it  is  more  convenient  to  treat  the  problems 
with  functions  of  complex  variables,  because  i  and  0  are  harmonic 
functions. 

First-order  approximation  of  the  canonical  form.-  It  has  been  known 
for  some  time  that  the  above  approximation  is  not  very  good  at  high- 
subsonic  velocity.  A  higher  order  of  approximation  can  be  made  as 
follows.  Instead  of  taking  the  approximation  as  given  in  equation  (31), 

the  approximation  can  be  made  accurate  to  the  slope  of  the  curve 
against  cd  at  the  free- stream  condition,  that  is, 


a  +  bcD 


(35a) 
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where 

co^)  (35b) 

dU3  /a^Doo 

b-( 

\  doo 

if  CD  =  0)^  is  chosen  as  the  free-stream  condition.  This  approximation 
is  shown  in  figure  10.  Then  equation  (19)  yields: 


^  13  djr  _  Q  (36b) 

9co2  3e2  a  +  to  8a) 


d20  820  +  b  M  -  0 

^2  S02  a  +  to  Scd  "  ' 


Introducing  cx>  =  p  +  ^  into  the  above  equations, 


-  — 

as?  +  ae2  ’  35  * 


(36b) 


(37a) 


S  +  ^+  iM  =  o  (37b) 

35^  ae2  “  ^ 


These  are  identical  to  the  equations  of*  Stokes  and  Beltrami  for  axially 
symmetrical  flow.  The  particular  solutions  of  interest  are  the  well- 
known  Bessel  and  trigonometric  functions,  and  are  not  given  here. 

These  differential  equations  are  invariant  under  translation  in  0. 
From  Lamb's  "Hydrodynamics,"  the  singular  solution  corresponding  to  a 

source  at  cjdq  =  -  or  o>0  =  0,  0Q  is 
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0(cd  +  £,e;  -  ±,e0)  = 


II"  *  W *  (•  -  »o)i 


1/2 


(38) 


Recently  Weinstein  (reference  56)  has  shown  that,  for  the  axially 
symmetrical  potential,  the  fundamental  solution  for  a  source  at 


(%  +  bBo)  is 


0(CD  +  |,0J  CDq  +  |,0O)  = 


ioge[(co  -  m0)2  +  (e  -  eQ)2]l/2 

*  f)1/2(»  + 1)1/2 


+  R 


(“ +  pe) 


(39) 


where  the  function  r(cd  +  is  regular  at  the  point  +  pS0) 

The  corresponding  singular  solution  to  the  stream  function  is 


6-6, 


(“  +  f)2  +  (6  “  e°)2 


1/2 


(40) 


The  fundamental  solution  of  the  stream  function  is  not  given  here, 
because  it  is  long  and  involved.  The  singular  solution  for  more  com¬ 
plicated  sources  can  he  built  up  easily.  Therefore,  if  the  boundary 
conditions  and  the  locations  and  types  of  sources  are  given  in  the 
hodograph  plane,  the  flow  in  the  hodograph  plane  can  be  found  with  the 
classical  technique  of  boundary- value  problems.  If  the  boundary  con¬ 
ditions  are  not  known,  the  series  solution  of  the  incompressible  flow 
may  be  considered  and  Chaplygin' s  procedure  followed. 

Second- order  approximat ions . -  Choose  the  following  approximation: 


where 


a  +  bco  +  cuf 


( 4la) 


2k 
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b  = 


ax1/2  q  a2K 1/2^ 


da) 


to2  /o>=a) 


(tec) 


c  = 


1  /d2K1//2 
2'  \  dco2 


(ted) 


CO=OX 


This  approximation  is  also  shown  in  figure  10.  Then  equation  (19)  yields 


a2t 

b  +  2ccd 

*1=0 

ae2 " 

a  +  1 xd  +  co£ 

ck D 

&  + 

b  +  2co) 

c 

II 

38 

ae2 

a  +  bo)  +  cocfi 

OO) 

(tea) 


(teb) 


These  equations  can  be 

cd'  =  - — - I— (1 0)  - 

(b2  _  4ac)1/2\ 


written  in  another  form  by  introducing 
-  j  ^if  b2  >  tec) . 


a2*  + 

aco’2 

+ 

am'2 


1  \<h 

CD1  +  1/aOD' 


=  0 


-3— )M-.0 

O)1  +  l/oO)f 


(tea) 


(teb) 


The  particular  solution  can  be  obtained  by  the  principle  of  separation 
of  variables,  if  it  is  assumed  that 
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+<0>',e)  «  )+2<v)(©)  -  (41(a) 

*<<»', e>  -  01(vV)02(v)<e)  -  sin  v!ls01(v)(co->  (44b) 

cos  vej  x 

where  and  0^  satisfy  the  ordinary  equations 


d2^ 

(dm')2 


d^ 

doo' 


(^5a) 


d2ii  2m'  d^l  2.  a 

- p  +  - 5 - —  "  v  ^  =  0 

(dm1)2  (m’)2  -  l  1 


(W 


which  are  special  cases  of  Stratton’s  equation  as  shown  in  reference  70. 
The  particular  series  solution  can  he  obtained  but  it  is  too  long  to 
give  here. 


Comments  on  above  three  approximations.-  If  the  above  three  approxi¬ 
mations  are  plotted  in  comparison  with  -^(log„  K1' 2 )  against  m  in 

dm '  e  ' 

figure  6,  they  are  not  satisfactory.  The  zero-order  approximation 
maintains  a  zero  value ,  although  the  true  value  of  -^^loge  K1/2  j  at 
may  not  be  zero.  The  first-order  approximation  becomes  roughly  the 


m. 


zero 


-order  approximation  to  -^-^loge  K^/2j  at 


m ,  while  the  second- 


order  approximation  becomes  roughly  the  first-order  one.  It  is  apparent 
that  the  above  procedure  can  be  repeated  to  find  all  the  approximations 
to  the  curve  in  figure  11,  but  this  kind  of  approximation  is  not  shown 
here.  Some  other  approximations  are  given  to  take  care  of  the  asymptotic 


behavior  of  —  [log  K1/2) 
doo  \  J 


as  go 


0  and  go 


00. 


2  6 


NACA  TN  2582 


Approximations  to  singular  "behavior  of  ^~(loge  Figure  11 

shows  the  curve  of  —flog  )  against  cd.  Some  simple  calcu- 

dco  \  e  ' 

lations  will  show  that  lim  ^^loge  K^2)  =  “  3"  an<^ 

co  • — >•  0 

A-jloge  ^(l  -  M2)1/2  — >0  as  co  - — >  00.  Therefore ,  a  simple,  hut 
good,  approximation  is  to  assume 


A(i„ge  W2)  =  ^ 


(46a) 


The  constant  a  can  he  chosen  in  two  ways: 

(l)  Choose  a  =  - -j  if  the  exact  behavior  at  sonic  velocity  is 
desired.  This  is  shown  in  figure  11. 


(2)  Choose  a  such  that 


,  if  the  flow  velocity  in 


the  neighborhood  of  free-stream  velocity  is  desired.  If  is 

infinitely  large,  it  automatically  reduces  to  Harman' s  approximation. 
This  case  is  also  shown  in  figure  11. 

This  approximation  is  important  because  the  fundamental  and  sin¬ 
gular  solutions  can  be  obtained  similarly  to  equations  (38)  and  (39) 

except  the  one-half  power  is  replaced  by  the  ^  power. 

The  next  approximation  along  this  line  is  to  assume  . 


(loge  K1/2) 


i(l  +  bu>) 


(46b) 


where  a  and  b  are  two  free  constants  to  be  chosen. 
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A  few  choices  can  he  made: 


l)  Take  a  = 


1  V-! 

—  and  b  to  be  such  that  — 


3  2K  “ooC1  +  toco) 

the  interest  is  in  the  sonic  and  free- stream  velocities. 


2}  Take  a  and  b  such  that  — — - - 

oo(l  +  boo) 


-(£), 


l— )  are  satisfied,  for  example.  This  case,  shown  in  the 

\2K/a^o.5mro 

graph  (fig.  11)  checks  well  with  the  exact  curve.  The  differential 
equations  are  simplified  to 


+<L± 


be2  <*(1  +  to)  ben 


(47a) 


_ 2_M,o 

^2  ^2  cd(1  +  bco)  8cn 


(47b) 


Similarly,  it  can  further  be  assumed  that 

£1  (46c: 

2K  oo(l  +  cod) 

where  a,  b,  and  c  are  free  constants.  One  case  is  shown  in  the 

graph  (fig.  ll)  with  a,  b,  and  c  determined  from  the  values  of  — 

2K 

at  0.500^,  oo^,  and  203^.  It  checks  very  well  with  the  exact  curve. 


3.2  -  Supersonic  Canonical  Forms  -  Equations  (23a)  and  (23b) 

Zero-order  approximation.-  Following  the  approach  of  Von  Karman  in 
the  subsonic  case,  there  can  be  chosen  in  the  supersonic  case  the 
approximation 


1/2  _  p, 


(M2  -  l) 


if/2 
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and 


Q  = 


la* 


(M2  -  l)^  ,  . 

- 5> - 34 


(k 8b) 


where  q  is  always  supersonic  in  the  domain  of  the  flow  pattern.  Then, 
the  coefficient  of  the  first  derivative  in  equation  (23)  is  zero,  and 


a2t 

bn2  be2 


(49a) 


0 

ha2  be2 


which  has  the  solutions 


(49b) 


*(n,e)  =  F+(fl  +  e)  +  G^(n  -  e) 


0.(n,e)  =  ^(n  +  0)  +  0^(2  -  0) 

which  are  well-known  classical  wave  functions.  The  approximation  is  a 
straight  line  parallel  to  the  ft- axis  as  shown  in  figure  12  with  ft  -  ft^. 

First-order  approximation. -  Along  the  same  approach  as  the  subsonic 
case,  it  is  easy  to  show 

yV  V 

bn2'  be2'  a'  +  Vft 


d20  d20  j  b’  M  =  0 

Sn2  8e2  a'  + 


.  (50a) 


(50b) 
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where 


(-K)1/2  -  0  a -£>l/g 


dft 


-Jft=ft 


b !  = 


a(-K)1/2 


dfi 


ft=ft. 


Introducing  Q  -  0,  +  —  into  the  above  equations , 


82t  _  8fi  _  1  =  0 

8n2  8@2  2  82 


(51a) 


820  820  1  80 

802  8e2  2  8ft 


(51b) 


If  the  following  characteristic  coordinates  are  introduced, 


\  —  ft  +  G 


[j.  =  -ft  +  6 


(52a) 

(52b) 


and  then, 


S^t  +  1  /St  St 

SA,  S|i  2(X  -  |j.)  \Sx  c )\i 


=  0 


(53a) 


820  1  /80  80 
8a.  8(i  2(x  -  n)  \8x  8|i 


=  0 


(53b) 
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which  are  Euler- Poisson' s  equations  with  P  =  P'  =  -  in  and  i 
in  The  Riemann  function  corresponding  to  \|r  can  he  shown  to  be 


X0,4q)  = 


(po  -  x)1//2(p  -  xo) 


1/2 


(h  -  X) 


i(-7’-§;H  (,k) 


where  Xq  (or  ftQ  + 
point  of  interest,  gF^ 


a 

and  |i0  correspond  to  ftQ  +  —  and  0Q,  the 
is  the  hypergeometrical  function,  and 


(X  -  X0)(n  -  nc) 
(x  -  M0)(n  -  >-0) 


It  becomes  too  cumbersome  to  write  the  Riemann  function  in  terms 
of  the  original  coordinates  ft  and  6.  This  approximation  is  shown  at 
ft  =  ftco  in  figure  12.  Since  the  Cauchy  data  are  difficult  to  assign  to 

the  potential  the  corresponding  Riemann  function  becomes  useless.  There¬ 
fore,  with  the  above  Riemann  function  $  the  stream  function  t  at  a 
point  (  ftQ,e )  can  be  found  by  the  Riemann  method,  if  the  Cauchy  data 

are  sufficiently  given.  In  general,  even  if  the  Cauchy  data  are  suf¬ 
ficiently  known  in  the  hodograph,  the  Riemann  method  is  rather  difficult 
to  apply. 


it 


Second-order  approximation. - 
is  not  difficult  to  show  that 


Following  the  same  procedure  as  before, 


82t  82i|r  (  1  +  1 

8ft2  8e2  \ft  +  1  &  -  1/dft 


+  I  _  1  +  _i—W  =  0 

8ft2  89 2  v>  +  1  "  l/dft 


(55a) 


* 


(55*>) 


5H 
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H. 


where 


H  = 


p— - — - — j-(a  -  21.)  (fb1  )2  >  Ite1  c') 

[(f>2  -  WcJ1/2'  2c';  V 


(56) 


2(-K)1/2 


dX2' 


(57) 


£2=fl„ 


The  Riemann  function  of  ^  is  under  investigation,  and  is  not  so 
easy  to  obtain  unless  a  long  series  solution  is  adopted.  However,  the 
particular  solutions  are  of  Stratton1 s  type,  but  are  not  given  here. 


Some  better  approximation  of  the  coefficient  —  log0  (-K)1/2 

dft  L  e  J 


can 


be  made,  but  the  solution  to  the  resulting  equation  will  be  more  diffi¬ 
cult  to  obtain. 


Comments  on  above  approximations.-  The  above  approximations  become 

worse  as  — ►  0  or  M  >  1.  Even  at  very  large  Mach  numbers,  they 

are  not  very  optimistic  approximations. 


3.3  -  Approximations  to  Transonic  Canonical  Forms  - 
Equations  (30a)  and  (30b) 

Take  the  transonic  canonical  forms  of  equations  (30a.)  and  (30b) 


8% 

c  s2* 

8e2 

cfgS  _ 

e  820 

ds2 

80  2 

1 

2 


M 


S0 


(30a) 


(30b) 
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Following  Guderley  (reference  30)  introduce  the  variable 


i  = 

(if-. 

_  6,2 

(supersonic) 

(58a) 

\2  /  e3 

ft2 

i  = 

(!)3^ 

_  e2 

(subsonic) 

(58b) 

W  63 

(JO2 

It  is  interesting  to  note  that 

£  >  0 

in  the 

supersonic  side  and 

t  <  0 

in  the  subsonic  side  of  the  hodograph  plane.  The  equation  £  =  1  corre¬ 
sponds  to  the  pair  of  characteristics  starting  at  6=0  on  the  sonic 
line.  The  independent  variable  9  can  be  eliminated  in  equation  (30a) 
and  equation  (30b).  Since  the  two  equations  are  equivalent  in  behavior, 
only  equation  (30a)  is  treated  here.  Thus  there  is  obtained 


Now  examine  what  conditions  must  be  imposed  on  e  — 

de 

so  that  the  variables  £  and  e  are  separable  in  the  above  equation. 
First,  assume 
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where  K  =  (^)2(l  -  M2)  is  a  known  function  of  e .  In  general,  the 

above  equation  cannot  he  solved  unless  the  following  conditions  are 
imposed: 


i-(loSe  (^)1  -  a 
3  31h  de  6  de  L.  e  \  e  /J 


s  T 


e  d*l  e  d[  (■ 
9  +i  de2  9  +3.  de  2  de  L  e  V 


=  a-aP, 


(62b) 


where  as  and  Ps  are  individual  constants  to  he  chosen  later.  Equa¬ 
tion  (6l)  can  then  he  written  as 


d2^o 

5(1  -  0  -^  + 


H  -  t1  +  a 


B  +  Ps)Sldf  "  asps^2  =  0  (63) 


which  is  a  hyper geometric  equation.  Its  singularities  are  £  =  0,  1, 

and  oo.  The  general  solution  about  £  =  0  is 


t) 


^2(0  = 


2Fllas  +  pPs  +  2} 


b  b  0? 


(U  I  <  1)  '  (64a) 


The  solution  about  £  =  1  is 


(1  -  02"as'Ps  2Fx(i  -  as,4  -  Ps;  §  -  as  -  ps  +  1; 


;  x  -  t) 


+2(0  = 


1/2  (1  -  5)2"as'Ps  2Fl(l  -  as,l  -  ps;  |  -  as  -  Ps  +  1;  1  -  5) 


(|l  -  i\  <  l)  (64 b) 
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Following  Kummer  (reference  6l),  all  the  24  solutions  about  the 
3  singular  points  £  =  0,  1,  and  »  can  be  obtained.  The  triple¬ 

valued  behavior  of  \|r2(0  for  £  =  0  can  be  shown  easily.  This 

agrees  with  the  results  of  Lighthill,  Guderley,  and  Carrier.  Under 
the  imposed  condition,  it  is  not  difficult  to  show  that 


Constant 


3(ps_tts) 

1  -  cse 


where  as  i  0S  and  cg  is  a  free  constant  to  be  chosen  later.  Also, 


if  by  definition  zs  =  ■-  •—  jloge  )j ,  its  approximate  value 
be  calculated  from  equations  (62a)  and  (62b): 


zg  can 


zs  =  1  -  3(|3S  -  <xs) 


1  +  cse 


1  -  cse 


3(Ps~as) 

3(3s-a.s) 


The  exact  value  of  zg  can  he  determined  from  equation  (25a)  with 

_ *  tr  _  /Po\2/,  x/2\ 


f(e)  =  €  and  K  =  (l  -  M2). 


It  is  necessary  now  to  choose  the  approximate  variable  z_  so  as 

to  have  the  same  ordinate  and  slope  as  zQ  at  6=0  (or  M  =  l)  if 

the  flow  at  sonic  speed  is  of  particular  interest.  It  can  he  shown 
that,  at  €  =  0, 


zs  =  0 


(67a) 


§£!>(# 


(67b) 
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These  conditions  determine 


2cse 
1  -  cse 


(68a) 


3S  =  as 


(68b) 


where  cs  sj  -O.525.  With  the  approximation  to  equation  (30a) ,  the 
approximate  differential  equation  becomes 


c)2i|r  2cs  cty 
8e2  be2  1  -  cse 

The  singular  solutions  of  this  equation  are 


(69) 


€ 


-  3a 
J  s 


(1 


*(e,e)  = 

€-3aS(i 


1.  9e£\ 

2’  4e3 j 


cse) 


-1  30 


2e 


3/2 


5.  3.  9S2' 


2  1 


6  2 


ke 


(70) 


where  as  should  be  chosen  to  avoid  the  limiting  line  in  the  physical 
flow  pattern.  It  should  be  noted  that  equation  (69)  and  its  equivalent 
equation'  of  0  (not  shown  here)  are  one  order  higher  in  approximation 
than  the  Tricomi  equation.  In  other  words,  the  Tricomi  equation  is 
equivalent  to  taking  zg  =  0  for  all  values  of  e.  For  negative 

values  of  as,  the  singularity  of  t  is  at  6  =  — .  If  <xs  is  posi- 

cs 

tive,  there  are  two  singularities  of  \|r,  one  at  €  =  0  and  other  at 
1 


6 
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4  -  DIFFERENT  APPROXIMATIONS  TO  CHAPLYGIN’S  SECOND 
DIFFERENTIAL  EQUATION  AND  THEIR  SOLUTIONS 


From  equations  (lOa) 


and  ( 10b ) , 


^  +  K  —  = 

b2  be2 


(71a) 


K  0  -  3?(l0ge  K)i  ’  °  (71b) 


is  an  implicit  function  of  a.  The  exact 

solutions  are  difficult  to  obtain.  Chaplygin  (reference  l)  in  his 
researches  on  subsonic  gas  jets  chose  K  =  1  so  that  his  equations 
become  Laplace  equations.  As  shown  in  figure  3  this  is  a  reasonable 
approximation  if  the  maximum  Mach  number  in  the  jet  is  much  less  than 
unity,  because  K  =  1  is  actually  the  asymptote  for  the  true  K  -  a 
curve  at  M  =  0.  Karman  (reference  7)  and  Tsien  (reference  52)  extend 
the  idea  to  the  linear  approximation  of  the  pres sure- volume  relation 
corresponding  to  the  free-stream  condition,  also  obtaining  a  Laplace 
equation  and  achieving  many  fruitful  results  of  technical  importance. 

In  the  transonic  flow,  of  course,  such  approximations  cannot  be  applied. 
Von  Karman  (reference  29)  in  19^7  gave  the  approximate  differential  equa 
tions  which  are  valid  when  the  flow  velocity  is  in  the  neighborhood  of 
sonic  velocity  and  when  the  body  in  the  flow  is  thin.  He  also  found  the 
transonic  similarity  law  which  gives  a  satisfactory  prediction  of  the 
wave  drag  for  thin  bodies.  In  the  hodograph  plane,  the  approximate 
differential  equations  are  of  the  Tricomi  type.  The  concept  of  the 
similarity  law  has  been  further  discussed  by  Kaplan  (reference  39)  and 
Guderley  (reference  31)  and  extended  by  Tsien  to  hypersonic  flow.  The 
main  contribution  of  the  transonic,  similarity  law  to  the  hodograph 
method  is:  Instead  of  investigating  the  flow  about  a  given  thin  body, 
the  flow  of  a  body  with  the  same  thickness  distribution  as  the  given 
one  but  of  vanishing  thickness  ratio  can  be  investigated.  Because  of 
the  vanishing  thickness  ratio,  the  boundary  conditions  can  be  simplified 
in  the  hodograph  plane  and  are  shown  to  be  consistent  with  the  approxi¬ 
mation  applied  to  the  differential  equations.  Under  these  conditions, 
the  problem  becomes  a  boundary- value  problem  of  the  type  studied  by 


ba 


where  K  =  ^  (l  -  M2) 
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Tricomi.  It  is  unnecessary  to  use  associated  boundary  conditions  from 
the  incompressible  flow  as  done  by  Chaplygin.  After  the  transonic  flow 
of  the  body  with  vanishing  thickness  ratio  is  obtained  the  transonic 
similarity  law  may  be  applied  to  find  the  aerodynamic  behavior  of  the 
given  thin  body. 

The  success  of  the  above  approximation  encourages  the  author  to 
seek  some  higher-order  approximations.  First  take  a  look  at  the 
K  -  a  curve  shown  in  figure  3.  It  has  two  asymptotes:  One  is  K  =  1 
and  the  other  is  cr  =  -0.2513.  It  is  a  monotonic  increasing  function 
of  a  within  the  range  -0.2513  <  a  <  at  cr  =  0,  K  =  0.  It  is 
understood  that  a  >  0  corresponds  to  M  <  1  and  a  <  0  to  M  >  1. 
The  main  interest  of  the  investigation  lies  in  transonic  flow.  There¬ 
fore,  if  possible,  the  approximation  should  be  so  chosen  as  to  maintain 
the  behavior  of  the  exact  K  -  cr  curve  at  and  near  the  sonic  velocity 
(a  =  0)  and  at  the  same  time  preserve  the  asymptotic  behavior  of  the 
exact  curve  as  M  — ►  0  (cr  — ►  »)  and  M  — *  »  (a  ►  -O.2513). 

Just  as  important,  the  approximation  equation  should  possess  solu¬ 
tions  within  the  reach  of  classical  mathematical  analysis.  With  these 
few  criterions  in  mind,  it  can  be  seen  that  the  Taylor  series  expansion 
of  K  about  a  =  0  is  not  a  favorable  choice,  although  the  first 
approximation  to  be  shown  below  is  of  this  nature. 


4.1  -  First-Order  Approximation  in  Neighborhood  of  Sonic  Speed 

Frankl  (reference  11)  in  19^5  obtained  the  Tricomi  equation  from 
equations  (71)  by  letting  K  2  K-j_  =  acr  and  consequently 


(72a) 


a20 

8a2 


+  acr 


820  1  80 
8e2  "  StT 


(72b) 


where  a 


and 


cr 


as  given  in  the 


earlier  definition.  The  comparison  of  the  approximation  with  the 

exact  value  of  K  is  shown  in  figure  3,  From  the  figure  the  range  of 


38 


NACA  TN  2532 


validity  of  this  approximation  is  very  narrow,  only  in  the  immediate 
neighborhood  of  the  sonic  velocity  a*.  ■  The  solutions  of  these  equa¬ 
tions  have  been  fully  investigated  by  Tricomi  in  his  famous  thesis 
(reference  53).  Carrier  and  Ehlers  (reference  57)  and  also  Weinstein 
(reference  40)  have  further  investigated  the  singular  solutions  of 
these  equations;  Tomotika  and  Tamada  (reference  33)  also  made  interesting 
contributions  to  the  singular  and  fundamental  solutions.  Recently 
Guderley  (reference  31)  made  an  extensive  study  of  the  singular  solutions 
and  showed  a  very  important  singular  solution  corresponding  to  a  family 
of  airfoils  at  sonic  speed.  Guderley  and  Yoshihara  (reference  36)  gave 
the  flow  over  a  wedge  airfoil  at  Mach  number  1.  They  employ  an  elegant 
method  of  attack  in  solving  the  problem. 

In  the  early  days,  Euler  and  Darboux  obtained  the  fundamental  solu¬ 
tions  of  the  same  equation  in  the  pure  elliptic  and  in  the  pure  hyper¬ 
bolic  domains  as  pointed  out  by  Tricomi  (reference  53).  The  importance 
of  Tricomi’ s  work  is  the  recognition  of  the  differential  equation  in  the 
mixed  domain,  that  is,  partly  elliptic  and  partly  hyperbolic,  and  the 
proof  of  the  existence  and  uniqueness  of  the  solutions  in  such  a  mixed 
domain.  Following  the  footsteps  of  Tricomi,  Frankl  (reference  11 )  has 
shown  the  existence  and  uniqueness  of  the  solution  of  Chaplygin' s  second 
differential  equation  in  the  mixed  domain,  particularly  on  two  problems, 
one  being  the  detached  shock  wave  of  a  wedge  and  the  other  being  the 
supersonic  jet  from  an  inclined-walled  vessel. 


4.2  -  Second-Order  Approximation 

From  the  comparison  of  the  first-order  approximation  and  the  exact 
curve  in  figure  3  it"  is  immediately  apparent  that  some  improved  approxi¬ 
mation  should  have  a  wider  range  of  validity  than  that  given  by  the 
linear  approximation.  The  usual  technique  of  taking  higher-order  terms 
in  the  Taylor  series  would  not  be  particularly  helpful  because  they 
would  not  improve  the  asymptotic  behavior  for  large  negative  values 
of  K  and  for  large  values  of  a.  Therefore,  choose 


K(cr)  =  K2(ct) 


acr 

1  +  c  a 


(73) 


where 


■  #)  ■  © 


any  one  of  three  ways: 


9.42  at  a  =  0  and  c  can  be  chosen  in 


6H 
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(1) 

(*■%' 

\  dcr2  ) 

1  W2/ 

(2) 

c  =  a 

so  that 

for 

transonic 

(3) 

c  =  -1/0.2513 

at  cr  =  0,  favorable  for  transonic  range 

co,  favorable 

— >  -0.2513, 


1  exactly  as 


that  K>>  =  K  — »  -00  exactly  as 
favorable  for  transonic  and  supersonic  range 


It  is  interesting  to  note  that,  in  equation  (73),  K? 


a 


00,  and  K> 


as 


-  — .  Thus ,  for  any 
c 


finite  value,  as  cr¬ 
one  of  the  three  choices,  K2  always  has  two  desirable  asymptotes. 

Therefore,  the  approximation  should  be  fairly  good  in  the  subsonic  range 
and  supersonic  range.  Of  course,  if  Chaplygin* s  procedure  of  using 
boundary  conditions  similar  to  those  for  an  incompressible  flow  is 
followed,  the  second  choice  is  a  favorable  one.  Cases  (l)  and  (2)  are 
shown  in  the  figure.  Both  seem  good  in  the  transonic  range.  Case  (3) 
is  not  shown.  The  only  known  second-order  approximation  is  given  by 
Loewner  (reference  73).  It  is  also  shown  in  figure  3  for  comparison. 


With  the  above  approximation  equations  (71a)  and  (71b)  become 


+  acr  _ 

5a2  1  +  CCT  c>02  " 

+  a a  +  1 

do2  1  +  co  ^02  <r(l  +  cor)  ba 

Equation  (74a)  can  be  solved  as  follows: 

Assume  the  variables  in  i|r(cr,0)  separable  and  let 


(74a) 


(74b) 


t(cr,e)  =  *i(n)(aHr2(n)(e) 


ko 
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where  the  positive  integer  n  is  introduced  to  characterize  the  nature 
of  the  solution 


*2(n)(e)  = 


c  2#  + 


(n  =  0) 


(nX£)  (n  ^  0) 

sm 


S  "I 

where  the  positive  constant  X  is  so  chosen  that  cog  (nk0)  is 

(n) 

periodic  for  any  fixed  interval  of  6 .  For  n  ^  0,  ^ ]_  (cr)  satisfies 

the  equation 


d2^(n)  .  t  =  0 

do2  1  +  Cd  1 


Introduce  a  new  independent  variable  z  =  2nX^j  ^tr  +  Obviously, 
in  the  range  -i  <  a  ^  «,  the  corresponding  range  for  z  is  6  ^  z  ^  „ 
Equation  (75)  yields 


d2yi(n)( 

dcr^ 


-*  'I*  lVl<n>(z)  *  0 


which  is  a  particular  case  of  the  Whittaker  equation  with  m2  =  i  (for 

l/2 

a  reason  to  he  shown  later).  Here  k  =  —  [—J  ^  0. 

2c  \c/ 


Before  the  choice  of  the  solution  is  made  from  the  known  results 
of  classical  mathematics,  a  clear  understanding  of  the  nature  of  the 
desired  solution  when  cr  — y  oo  or  when  the  flow  becomes  incompressible 

is  necessary.  As  cr  — >  00  >  z  — > 00  and  \|r^  ^(z)  must  behave 
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kl 


+4 


like  e  ^  from  the  limiting  case  of  equation  (76).  Furthermore 
sin  n X6 

will  be  just  a  single  term  of  the  convergent  infinite 

(n) 


cos  n\0 

series  which  represents  to  be  shown  later.  The  term  ijr^n>,(z) 

should  remain  finite  when  a  — >  00.  Therefore,  \|f^(n)(z)  ~  e  2  as 
a  — ^  00  is  necessary. 

With  the  desirable  asymptotic  behavior  of  ty-^n\z)  at  large 
values  of  a  in  mind,  the  only  choice  is 


=  wK  (z) 
k’~2 

m 

whose  integral  representation  (see  reference  6l)  is 


W<n>  (z) 
k  -A 
2 


r(k  +  l)e  2  zk 


2«i 


P(Of) 


(-t) 


-k- 


H1  *  r 


e_t  dt 


(77) 


where  the  path  of  integration  is  a  contour  in  the  complex  t-plane 
starting  from  00  just  above  the  real  axis,  encircling  the  origin  in 
the  positive  direction  (counterclockwise),  and  returning  to  the  starting 
point  just  below  the  real  axis.  In  general,  for  k  >  0  there  are  two 
branch  points,  one  at  t  =  0  and  the  other  at  t  =  -z,  a  point  on  the 
negative  real  axis  (since  z  >  0).  Thus  the  path  of  integration  must 
be  chosen  so  that  it  will  not  encircle  the  branch  point  t  =  -z. 

The  asymptotic  expansion  for  large  values  of  a  (not  simultaneously 
of  large  values  of  n)  is 


W(n)!  (z) 
k’'I 


2  zk 


1  + 


OO 

L 

S=1 


-  ■  *  d! 


~  e 


s!zs 


(78) 
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For  the  case  of  large  values  of  n,  the  asymptotic  expansion  of 
W(n)1  (z)  is  much  more  complicated  owing  to  the  fact  that  both  k 

k » —  g 

and  z ,  being  proportional  to  cr,  become  large  simultaneously.  Very 
little  information  on  such  expansions  is  available. 

For  numerical  calculation,  the  corresponding  series  solution  is 
desirable.  Jeffreys  (reference  71 )  condenses  early  developments  of 
Whittaker,  Goldstein,  and  Stoneley  and  puts  into  compact  form  the 
relation  between  the  Whittaker  function  and  the  confluent  hypergeometric 
function.  The  following  discussion  closely  parallels  his  work  except 
for  some  changes  in  notation.  For  2m  =  ip  (p  is  a  positive  integer), 
the  Whittaker  function  can  be  represented  by  a  combination  of  Kummer's 
series  when  the  limiting  value  is  taken.  Thus, 


Wk  +u  (z)  =  lim 

2m — >±p 


r ( -2m)  .  .  T(2m) 

M,  (z)  +  — - - M.  (z) 

Tc,m  /I  ,.\  K,-m 


(§-“-*) 


(i  +  «  -  k) 


where 


-Z  +m  +  i 


Mk,±m  (z)  =  6 


_  ^  2  ~~  2  r(±2m  +  1) 


“  ri+m  -  k  +  i  +  s)  zs 


(±»  -  k  ♦  |)  ^  r<±2“  +  1  +  s)  s: 


=  e 


±m  + 


ifi( 


±m  -  k 


+  i  ,±2m  + 

2.’ 


l; 


2®\(a,  y, 
=  m  - 


z)  being  the  confluent  hypergeometric  function, 
k  +  i,  7  =  -2m  +  1  according  to  Jeffreys'  notation- 


~  Z  m+  1  V 

W  (z)  =  lim 

e  2  z  2  U|-m  -  k  +  -,-2m  +  1;  z) 

k,+p  2m->+p 

L  \  2  ;J 
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where 


U 


(-m  -  k  +  i,-2m  +  lj  zj  =  ^  ^ - -  ^(m  -  k  +  i,2m  +  1;  z) 

r  (-m  -  k  +  JJ  x  ' 


r(2m)z2m  „  /  ,  1  \ 

7 - IY  iFi("m  -  k  +  p-2*  +  i;  z) 

r(m  -  k  +  |J  '  2 


He  gives  the  solution  of  the  limiting  case  2m - >-  ±|u.  For  the  present 

case  if  m  is  chosen  equal  to  -i  (7  =  2  in  his  notation),  there  can 
he  written 


W  ]_  ( z)  =  e 

k,-? 


zju^z)  +  U2(zj] 


(79a) 


where 


u  (z)  =  - i -  z"1 

r(i  -  k) 


(79-b) 


u2^z)  =  l1^1  ■  k >2’’  z)  (^°Se  z-  FU)  -  F(°)  +  P(-kj]  +  u3  (79c) 


The  function  F(£  )  is  called  the  digamma  function  and  is  generally- 
represented  by  \|r(£  +  l)  (reference  74). 

F(0  =  ^0ge  r(^  +  13 


lim  logp  n  - - - - 

n — >00  V  £  +  1  £  +  2  5  +  3 


(80) 


(F(0)  =  Euler  constant  (-*7)  =  -0.57722.) 


u 
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For  the  present  case, 


r(-k) 


(l  -  k)  (2  -  k)  .  ■  .  (s  -  k) 
s’(s  +  1): 


r  -  k 


2VI*  1.-3- 

/ .  r  ■  s  +  1 

r=l 


”  .  x  /  s  s 

-Y  rd  <-  3  -  k)  .b  /y  _i —  9  1 ...  1  — 3 

frj  r(-k)r(l  -  k)  r  -  k  8  +  1 


(k  ^  integer)  (8l) 


Then,  the  series  solution  of  equation  (76)  is 


00 

i(or,e)  =  (c^cr  +  c2)(c30  +  04)  +  AnVT »  1  (z)  sin  (n\0  +  c^)  (82) 


where  c,  c^,  c^?  c,  ,  A  ,  and  a  are  constants  to  he  determined 

V  2  3  n  n 

from  the  given  boundary  conditions  or  from  the  boundary  conditions 
associated  with  an  incompressible  flow  as  in  Chaplygins  method,  if 
necessary.  Of  course,  the  convergence  of  the  series  must  be  established 
for  the  particular  problem  in  order  to'  be  sure  the  above  representation 
is  correct. 


4.3  -  Third-Order  Approximation 

The  K  -  a  relation  can  be  approximated  more  closely  by  assuming 
an  analytic  function  of  a  in  the  range  -  i  <  cr  <  00  of  the  form 


K  =  Ko  = 


ac(l  +  bcr) 
(1  +  ca)2 


(83) 
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such  that 

(1)  K3  =  K  =  0  at  g  =  0 


(2)  S  =  ffl  .  a  at  „  .  0 
da  da 


,  x  d2Ko  d2K 
(3)  — if-  =  — -  at  a  =  0 
d ad  d ad 


(4)  K^( cr)  also  possesses  two  asymptotes,  one  of  which  can  be 

made  to  coincide  with  an  exact  asymptote  of  K(g). 

The  conditions  (3)  and  (4)  determine  the  constants  b  and  c.  The 


above  approximation  with  Ko  — >  —  =1  as  a 

3  2 


oo  is  shown  in  fig¬ 


ure  13.  It  checks  very  well  with  the  exact  value  of  K  for  the  range 
from  subsonic  through  transonic  up  to  supersonic  regimes.  The  other 
asymptote,  a  =  -O.2583  as  £3  =  does  not  differ  greatly  from 

a  =  -0.2513  for  the  exact  value  of  K.  There  is  another  important 
advantage  of  this  choice,  because  the  boundary  conditions  of  incompres¬ 
sible  flows  can  be  borrowed  as  in  the  Chaplygin  procedure. 

Introducing  equation  (83)  into  equations  (71a)  and  (71b), 


crt  +  ag(  1  +  bg)  cT^  = 
bo2  (l  +  eg)2  80^ 


(84a) 


and 


820  +  ag(l  +  bg)  820  (2b  -  c)g  +  1  80 

bo2  (l  +  eg)2  80  2  or(l  +  bg)(l  +  eg)  8g 


(84b) 


To  solve  equation  (84a),  assume  that  the  variables  are  separable  and 


*(n)(g,e)  =  ^1(n)(g)i2(n)(e) 
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where  the  superscript  (n)  is  used  to  show  each  solution  is  related 
to  n,  a  positive  integer, 


Vn'(e)  = 


(nxe)  (n  ±  0) 


c30  +  c^  (n  =  0) 


and  ^  n  (cr)  satisfies 


d^l^  \  _  n2X2aq(l  +  bg)  ^  (n) 
dg2  (l  +  eg)2 


With  the  introduction  of  a  new  independent  variable  z  =  2nA  M(i  +  ccr), 

c2 

equation  (86)  can  he  transformed  to  the  well-known  Whittaker  equation 


ia2  V  '  z2  / 


where 


k-  0 

2c2  Vi 


(c  >  2b 


(88a) 


1  n2\2a(c  -  b 


imaginary  if 


i2A.2a(c  -  b)  .  1 


i2X.2a(c  -  b)  1 


(88b) 


real  if 


7H 
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The  fundamental  system  of  solutions  of  equation  (86)  is 

(n)  „  ,  (n)  .  .  .  , 

cnWk,m  (z>  *  (-*)  <”  +  °> 

c-j^cr  +  c2  (n  =  0) 


ti(n>(z)  -  \ 


(89) 


where  Cn  and  Dn  are  constants  to  he  determined. 

The  Whittaker  functions  can  he  expressed  in  terms  of  the  confluent 
hypergeometric  functions  and  Rummer*  s  functions  if  m  /  ±  ■—  as  follows 


w[n)  (z) 
k,m 


F(-2m)  (n)  ,  v 

—  M.  (z)  + 
l\  k,m 


(-"  - k  + §) 


r<2B,)  <2>  (90a> 


r(m  +  k  +  I) 


and 


«<?!*  <-■>  •  ,  r(~2m)  ^  (-«)  *  <-> 


^-m  +  k  +  ij  '  r^m  -  k  + 

(|arg(-z)|  <  |  jtj 


(90b) 


where  the  Kummer’s  series  are 


—  ±m+i  t* f j. o™  .  u  Jg_  r^tm  -  k  +  g-  +  sj 


„(n)  f  *  "2  -m'r  2  r(±2m  +  l)  V-  '  V - “  '  Z  '  “/ 

Ml  '  (z)  =  e  z  — - ^  \  i - 

’  m  p'±m  -  k  +  4 1  TZr\  r(+2m  +  1  +  s)s! 


-f  ±m+i 

=  e  z  ^ 


(±m  -  k  +  I)  r(i2m 

lFl(±m  -  k  +  p±2m  +  1;  zj 


M^(-z)=eiz±m+i  r(±2m  +1) 


■k,±m 


(±m  +  k  +  £  +  s) 


z) 


+m+  i 


=  e^  z 


r^+m  +  k  +  r(±2m  +  1  +  s) 

■jF^tm  +  k  +  i,+2m  +  1;  -zj 
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Equation  (92a)  gives  an  exponentially  damping  function  of  z  and  equa¬ 
tion  (92b)  is  an  exponentially  increasing  function  of  z  which  is 
divergent  as  cr  — Therefore,  the  only  choice  is  Dn  =  0.  Thus, 
there  can  be  written 


00  r~  .  . 

iKct,0)  =  (cjff  +  02)^036  +  c]j  +  AnWk^m  (z)  cos  nX0  + 


B  (z)  sin  nX6 

n  k,m  v  '  J 


where  cn ,  c~,  c_,  c.,  A  ,  B  ,  and  X  are  determined  from  the 

1'  2  3  A7  n7  n7 

given  boundary  conditions  or  from  the  boundary  conditions  of  incompres¬ 
sible  flow,  if  necessary.  If  An  cos  and  Bn  sin  0^  are  introduced,, 

there  can  be  written 


00 

i(o,0)  =  (cxa  +  c2^c3e  +  clj  +  ^  AnWk,m  (z)  sln  (nX0  +  °o) 


It  is  interesting  to  note  that,  when  m  =  +i,  equation  (87)  reduces  to 

equation  (76).  Thus  the  second-order  approximation  is  just  one  par¬ 
ticular  case  of  the  third-order  approximation  .with  b  =  c.  Of  course, 
when  the  choice  is  made  that  b  =  c  =  0,  both  cases  reduce  to  the  first- 
order  approximation,  equation  (72a). 

Along  this  line  of  thought,  if  there  is  imposed  k  =  0  or  2b  =  c, 
equation  (82)  becomes 


K  = 


The  differential  equation  becomes  the  normal  form  of  the  Bessel  equation 


2,ir  (n) 


1  __2 
T  -  m 

1  +  it _ 

4  *2 


(n)  = 
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where 


-  ^  H?(i + -) 


(97a) 


m  = 


1  n^X^a 


2c  ^ 


Its  solution  can  be  written  as 


(97b) 


♦ 


i(n)C*)  ->££(*)  ■  |f «x»(n) |) 


where 


is  the  modified  Bessel  function  of  the  second  kind  with 


complex  order  m.  The  function  and  no  others  satisfy  the 

requirement  of  boundedness  in  value  as  0  — y  00  •  This  can  be  shown 
by  the  asymptotic  expansion 


(98) 


There  is  another  interesting  feature  of  this  approximation  if  the 

cLK 

following  conditions  are  chosen:  a  =  —  at  a  — ►  0  and  c  such  that 

da 

c  =  —  or  K  — >  1  as  a  — >  00 .  Then, 

2 


(99) 
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and 


( 100a) 


m 


1  4n2\2 

*  "  a2 


(100b) 


The  advantage  of  this  is  that  the  boundary  conditions  can  still  he 
brought  from  the  incompressible  flow. 

The  accuracy  and  usefulness  of  this  third-order  approximation  is 
discussed  in  the  last  section  of  the  paper. 


5  -  HYPOTHETICAL  GAS  LAW  CORRES PONDING  TO  APPROXIMATIONS 
OF  CHAPLYGIN *  S  SECOND  EQUATION 


In  equation  (71a),  if  \|r  =  ^(cr)^^)  is  introduced;  the  principle 
of  separation  of  variables  yields: 


a2ti 

da2 


n2X2K(a)i|r1  =  0 


(101) 


where 


K  = 


1  -  d]7 

(l-T^l 


(102) 


a  =  -0.2513 


(1  -  t) 


1/2 


1  - 


1  -  1  +  (1  -  t)2 

3  5 


+  tanh"1  (1  -  t)1/2 
(103) 


d£ 

dq 


_P_  1 
Po  <3 


(104) 
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Since  K(cr)  is  an  implicit  function  of  a,  some  approximations  of 
different  orders  have  been  made  in  the  preceding  section  to  obtain 
useful  solutions  of  tp.  It  is  natural  to  ask  what  kind  of  hypo¬ 
thetical  gas  laws  will  correspond  to  the  different  approximations. 

It  is  understood  that  a  is  just  another  way  of  expressing  the 
velocity  magnitude;  consequently  it  is  defined  once  and  for  all  time 
by  the  exact  gas  behavior  as  shown  in  equation  (103).  The  variable  a 
is  independent  of  the  hypothetical  gas  law  for  each  approximation 
shown  in  the  preceding  section.  The  problem  now  is  to  find  the  func¬ 
tional  relations  of'  p  and  p/pQ  in  terms  of  cr  as  the  independent 
variable.  First  of  all,  by  definition, 


^  =  (a*)2  (105) 

dp 

and  the  differential  form  of  Bernoulli’s  equation  gives 


dp 

dq 


=  -pq 


If  equation  (104)  is  rewritten 


(106) 


it  can  easily  be  shown  that 


M2  = 


d,/p 
da 


(107) 


Consequently,  if  is  expressed  in  terms  of  q_1  and  its  derivatives 

with  respect  to  a,  equations  (106)  and  (107)  yield 


M2  =  1 


2 


(I08a) 
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On  the  other  hand,  equations  (102)  and  (106)  give  an  alternate  expres¬ 
sion  for  M2  in  terms  of  q”-*-  and  its  derivatives 


M2  =  1 


(S 


K(ct)  =  1  - 


^(l06e  -a"1) 


K(ff) 


( 108h ) 


Equating  the  above  two  equations,  a  differential  equation  is 
obtained  of  q-1  with  respect  to  a, 


dg(q-l) 

da2 


K(a)q--*-  =  0 


(109) 


which  is  the  same  equation  as  equation  (101),  if  q-1  is  placed  for  \|r 
and  nX  is  set  equal  to  1.  If  the  solutions  of  equation  (101)  are 
known,  the  solutions  of  q  can  similarly  be  obtained  explicitly  in 
terms  of  a. 

If  now  the  solution  of  equation  (109)  is  substituted  into  equa¬ 
tion  (106);  the  density  ratio  is  obtained  as  a  function  of  a.  Dif¬ 
ferentiating  this  expression  with  respect  to  a  as  given  in 

a.  The 
has  to  be 


(no) 


the  relations 
those  already 


K(cr)  in  order 
approximate 

solution  of  q"1  in  terms  of  a  called  q_1(a).  With  the  term  q_1(a), 
p(cr)/p0  can  be  obtained  from  equation  (104)  and  p(a)  from  equa¬ 
tion  (110).  Then,  the  approximations  can  be  compared  with  the  exact 
values,  both  in  terms  of  a . 


equation  (107),  the  Mach  number  is  obtained  in  terms  of 
derivatives  above  involve  only  differentiation  but  p(a) 
obtained  by  the  integration 


=  /p< 


,q2  da 


Of  course,  if  equation  (109)  can  be  solved  exactly, 
of  p/pQ  and  p  with  respect  to  a  will  coincide  with 
obtained  from  the  exact  gas  law. 

Now,  introduce  different  orders  of  approximation  to 
to  solve  equation  (109).  Then,  obtain  the  corresponding 
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In  the  preceding  section,  the  best  approximation  to  K(cr)  is  the 
third  approximation 


acr(l  +  bcx) 
(l  +  ccr)^ 


(111) 


*■ 


where  a,  b,  and  c  are  constants  already  chosen  according  to  certain 
considerations.  Substituting  into  equation  (109),  there  results  the 

differential  equation  of  q-1  corresponding  to  the  approximation 
of  q"1, 

_  ag(l  +  boQ(--l)  =  0  (112) 

dcr^  (l  +  ccx)^ 


The  proper  choice  of  the  solution  is 


q-1 


EM.  (z) 
k,mv 


+  FM, 


(z) 


(113) 


where 


(l  +  ccr)  >  0 


>  (11*0 


m 


The  terms  Mv  (z)  and  Mv  ^(z)  are  Kummer's  confluent  hypergeometric 
ic  ^  mx  )  — m 

series  (reference  6l)  as  shown  in  equations  (90a)  and  (90h)  and  E 
and  F  are  arbitrary  constants  to  be  determined  from  the  following 
boundary  c  ondit ions : 


8H 
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(a)  At  a  =  0,  q  =  q  =  a* 

(b)  At  a  =  0,  1 


d(q-!) 

dq-l  _ 

te>\  . 

1 1 

(7  +  1' 

)p 

dor 

da 

\pqLo 

a*' 

(  2  / 

) 

must  go  to 

00  as 

a  — y  00 

(M 

q  — 

0). 

The  above  are  the  correct  boundary  conditions  because  it  is 
intended  to  choose  the  hypothetical  gas  law  to  agree  with  the  exact 
one  at  sonic  velocity  (cr  =  0)  to  the  highest  possible  order  of  approxi¬ 
mation  and  to  preserve  the  asymptotic  behavior  as  a  — y  oo. 

2  y 

If  there  is  introduced  z  =  zQ  =  — ' —  at  a  =  0, 

c2 


a*zQW 


(z0)  '  C  £  ) 


-  M,  (z0) 
c  k,-m\  °J 


(115) 


a*zQW 


-  *cK,Jzc) 


(11 6) 


where  v(*0)  -  ^(soK^o)  ’  V*M  W  is  the  Wronskian. 
From  equation  (104),  the  approximate  density  ratio  is 


1  1  EMk,m(z)  +  FMk,-m(z) 

Po  ‘  zQc  El&>a(z)  +  FM^_m(z) 


(117) 


which  is  plotted  in  figure  l4.  It  differs  very  little  from  the  exact 
curve.  From  equations  (108a)  and  (108b),  the  approximate  Mach  number  is 


M  =  1  - 


(118) 
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which  is  shown  to  compare  reasonably  well  with  the  exact  curve  in 
figure  15. 


Since  the  hypothetical  gas  law  differs  so  little  from  the  exact 
value,  it  seems  justifiable  to  use  the  exact  gas  law  to  replace  the 
hypothetical  gas  law  if  necessary,  particularly  in  the  neighborhood 
of  the  sonic  velocity. 

To  the  author's  knowledge,  the  only  available  high-order  approxi¬ 
mation  is  Loewners'  approximation  (reference  73).  It  is  the  basis  of 
Carrier  and  Ehlers’  investigation  on  channel  flow  (reference  32).  For 
comparison,  his  approximation  in  gas  behavior  is  given  in  figures  l4 
and  15.  His  approximation  is  correct  at  the  sonic  velocity  to  the 
second  derivative  of  the  K  -  a  curve. 

For  the  present  second  approximation  b  =  c  or  m  =  -  is 

chosen.  When  2m  is  an  integer  only  one  of  the  series  solutions  is 
valid,  namely,  Mkj_m(2)-  If  the  solution  Mk^_m(z)  is  retained  then 

the  second  independent  solution  can  be  obtained  by  using  the  limiting 
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value  of  Mj^_m(z)  as  m  - i.  It  is  quite  similar  to  the  develop- 

ment  of  the  preceding  section.  No  further  details  are  necessary. 

For  the  case  2b  =  c,  k  =  0.  Similar  solutions  can  be  obtained 
with  slightly  different  boundary  conditions. 

6  -  TRANSFORMATION  BETWEEN  HODOGRAPH  PLANE  AND  PHYSICAL  PLANE 


With  the  introduction  of  dcr  = - -  —  the  differential  equations 

P  q 

for  and  0  from  equations  (10)  and  (l 6)  become 


M  =  K(tr)  ^ 
8a  80 


( 121a ) 


M 

de 


ae 


(121b) 


where  K(cr) 
in  \|r  is 


As  shown  before ,  the  differential  equation 


—  +  K(  a)  =  0  (122) 

8a2  8e2 


with  the  exact  value  of 


K(cr)  approximated  by 


K3(a)  = 


aa(l  +  ba) 
(l  +  ca)2 


Introducing  t(a,0)  =  ^(a)^©). 


(123) 
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where  n  is  an  integer  and  X  is  positive  and  real;  A,  is  introduced 
here  in  order  to  make  ^2(nXe)  =  \|r2(nX0  +  2n)  periodic.  For  simplicity 

X  =  l  is  chosen  in  the  later  treatment.  With  the  introduction  of 


2n 


fat>(l  +  ccr)  (not  to  be  confused  with 


presently) ,  the  general  solution  of  equation 
be 


Z  to  be  introduced 
(122)  has  been  found  to 


► 


r 


t(d,e)  =  A0(cr  +  ff0)(e  +  e0)  +  ^  Vk?!  (z)  cos  ne  + 

n=l 


y  Bnwk^m  sin  n6  (12^) 

n=l 


where  AQ,  An,  Bn,  aQ,  and  Qq  are  to  be  evaluated  from  the  given 

boundary  conditions  for  a  particular  problem  at  hand,  or  from  the 
solution  of  the  corresponding  incompressible  flow.  The  superscript  (n) 
is  applied  to  the  Whittaker1 s  function  in  order  to  show  its  relation 
to  n.  Suppose  that  all  these  constants  are  known  and  the  right-hand 
side  of  equation  (124)  is  assumed  convergent  and  represents  y(cr,0); 
it  remains  to  transform  the  results  obtained  in  the  hodograph  plane  to 
the  physical  plane  so  that  the  problem  is  solved  in  the  physical  plane. 
The  procedure  of  carrying  out  such  a  transformation  is  given  here. 

First  introduce  q  and  0,  the  inclination  of  the  velocity 
vector,  as 


u  =  q  cos  0 
v  =  q  sin  0 


50  _  po  5^ 


where  u  =  —  - - is  the  local 

5x  P  5y 


and  v  =  ^  ^  is  the  local 

5y  p  5x 


.  50  _ 

-  _  -  — -  -UU  UiiU  -L.'-'  U 

p  8x 

Then  the  total  differentials  of  0, 


velocity  component  along  the  X-axis 

velocity  component  along  the  Y-axis, 
ijr  are 
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d0  =  M  dX  +  ^  dY 

8x  8y 


=  q(cos  0  dX  +  sin  0  dY) 


(125a) 


d*  =  M.  ax  *  M  dY 
8x  8y 


=  —  q(-sin  6  dX  +  cos  6  dY) 
Pq 


(125b) 


Or,  in  short. 


dZ 


(125c) 


where  Z  =  X  +  iY  is  a  complex  variable  in  the  physical  plane.  With 
cr,  9  as  independent  variables  in  the  physical  plane, 


d 0 


M  da  + 


(126a) 


di 


**  dc  ♦  M  de 


8a 


80 


(126b) 


The  derivatives  of  0  can  be  eliminated  entirely  if  the  relations 
given  in  equations  (121a)  and  (121b)  are  introduced  in  equation  (126a) 


d0  =  K(c)  p-do  d 0 

80  8a 


(127) 
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Substituting  equations  (l26b)  and  (127)  into  equation  (l25c), 
equation  (125c)  can  be  rewritten  as  the  exact  differential 


,  ^  8z  _  8z  * 

dZ  =  —  dcr  +  —  d0 

8a  80 


or 


da  + 


8^  .  po  8^ 

8a  p  80 


(128a) 


8z 

8a 


8t 

80 


i  £0  8t 

p  8a_ 


(128b) 


8z  _  e±e  /  8t  +  i  po  8^  \ 
80  q  \  8a  P  80  / 


(128c) 


The  above  relations  are  conformal  and  uniquely  exist  as  long  as  the 
Jacobian  determinate 


0  or  .  (129.) 

8(o,e) 


In  the  case  of  Y )  _  q  (x,Y)  some  singularity  may  occur  in 

8(a,0) 

the  physical  plane.  If  —  =  00  some  singularity  may  occur  in  the 

8(a;0) 

hodograph  plane  as  shown  in  reference  48. 


Before  carrying  out  the  calculation,  a  few  symbols  are  introduced 
to  make  the  presentation  a  little  clearer.  Let 


I(r(0)  =  {a  +  aQ)(e  +  0O) 

(z)  cos  n0  =  cos  n@ 

t  (n)  =  w[n)  (z)  sin  ne  =  w(n) 
s  K,m 


(130) 


sin  n0 
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Also  represent  Z  =  Z^0)  corresponding  to  a: 

\|^n^(z).  For  future  use,  rewrite  equation  (106)  as 


and  Z 


=  7.(n) 


Po 
q  — 
p 


(131a) 


i 

p  da  p 


(131b) 


±(?o 

dcrlpq 


(131c) 


To  carry  out  the  evaluation  of  Z  in  terms  of  a  and  0  is 
rather  too  long  for  the  space  available  in  this  report.  The  basic 
method  of  calculation  is  rather  simple.  Since  in  equation  (128a) 

dZ  =  7— -  da  +  d©  is  an  exact  differential,  t—  (or  t—)  can  be 
5a  50  9  50  V  5a/ 

integrated  to  obtain  Z  except  for  an  unknown  function  F(a)  (or 
F(0))  to  be  determined.  The  unknown  function  F(a)  (or  F(0))  can 


be  evaluated  from  the  remaining  relation 


as  (or  as') 

5a  \  50/ 


Of  course,  in 


the  operation,  repeated  use  has  been  made  of  the  differential  equation 

Q 

d  ( a )  p  .  . 

- — r - n  Ktp(cr)  =  0  or  Whittaker*  s  equation.  For  further  details 

da2  * 

the  following  section  may  be  consulted. 
n=  0. -  In  this  case, 


z(o)  .  M  t  1T(o) 


where  Z, 


^  -  y(°  +  ffo)  -  x]  -  i(e  +  eo)  -  +  zc 


(o)  =  Y  (o) 


xow;  +  iYo 


is  the  constant  of  integration. 


(132) 


MCA  TN  2532 


n  =  1.-  For  n  =  1,  the  function  F^^(cx)  is  an  integral  as 

F(1)(0)  ,i  Ilk  5  W<1>1  d, 

d  ,  pq  dcr 


(133) 


The  term  Z^-1'  can  expressed  by 


z  (1)  _  z  (1)  -  ie2ie  dW*1)  h  P_c 
C  C°  4q  L  dcr  P 


6 

dW^1) 

2q 

—  da 

+ 


*  -1  pn  dW'--L'  (l) 

i  i  12 - K VIK  !  do 


2  q  p  do 


z  (1)  _  Z  (1)  =  £!^_+  Po  w(i)  _  1£  dW^_  £0  w<i)  + 
s  so  4q  da  p  2q  da  p 


1  1  Po  dW1 


KWV  da 


2.  qlp  <*a 


(I3to) 


where  Z , 


and  Zs  are  the  integration  constants. 


o 

n  /  1  or  0,-  For  the  cases  n  /  1  or  0,  the  function  F^n^(a) 
found  to  be  a  constant 


Z  (n)  _  z  (n)  =  ie1(n+1^  dW^n)  +  ^o  w(n 
c  C°  2(n  +  l)q  do  p 


ie-i(n-l)8  dW(n)  ^  n£o 
2(n  -  l)q  _  do  P 


( 135a) 


9H 
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In  more  detail,  the  expression  is 


Z  -  ZQ  =  — 


♦  •„)  - 1  -  he  *  O]  *  (lAl 


hi  „U)  ( J  .  5>  w*1*  (z) 


1A-,  +  B-,  .  P  (z)  (1  \ 

— - -  -  -2.  — ^ - K(a)wi1^  (z)  da  + 

2  q  p  dcr 


K  +  Bn)'1<n+l)eK“i  <*>  np0  (n) 


2(n  +  l) q 


—  WH  (z)  + 

P  k'm  J 


(-lAn  +  Bn 


+  Bn)e-i(n-l)0pi“i  (2)  nP0  u(n) 

2(n  -  l)q  da  P  k’E 


Wi  '  ( z ) 

k,m  '  ' 


(137b) 


2n 

where  z  ~  ~  ^*^(1  +  ccr)* 

The  above  result  is  also  true  for  the  cases  of  the  second  approxi¬ 
mation  h  =  c  /m  =  -  or  of  the  approximation  h  =  2c  (k  =  0). 


7  -  FLOW  OF  COMPRESSIBLE  FLUID  THROUGH  AH  APERTURE  OF  A  TWO- 
DIMENSIONAL  INCLINED-WALLED,  STRAIGHT-EDGE D  NOZZLE 


To  apply  a  critical  test  to  the  present  investigation  it  seems 
reasonable  to  compare  the  present  approximation  with  a  well-known  flow 
pattern  which  has  been  studied  by  early  explorers  using  the  exact  method. 
For  instance,  Chaplygin  gave  an  application  of  his  investigation  to  the 
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efflux  of  a  gas  from  an  infinite  vessel.  Later  Lighthill;  following  his 
approach^  repeated  the  example.  Both  considered  the  particular  case 
with  the  walls  normal  to  the  jet  and  both  brought  the  boundary  condi¬ 
tions  from  the  incompressible  flow  and  treated  the  maximum  jet  velocity 
up  to  the  sound  velocity. 


The  present  application  gives  a  more  general  type  of  such  flow  - 
flow  through  an  aperture  of  an  inclined- walled;  straight- edged  nozzle. 
Besides;  the  problem  is  treated  directly  as  a  boundary- value  problem. 


The  question  of  the  maximum  velocity  in  the  jet  depends  on  the 
value  of  the  ratio  of  the  pressure  p 


surrounding  the  jet  to  the  stag- 

7 


p  (y  +  l\  7“1 

nation  pressure  pQ  in  the  vessel.  As  long  as  —  ^  \  ~o — )  >  the 

P0  \  d  / 

maximum  velocity  can  never  exceed  the  sound  velocity.  The  boundary 
value  of  the  stream  function  ^  is  clearly  defined;  therefore  it  is  a 

7 

7=T 


direct  boundary- value  problem.  For  the  case 


the 


velocity  of  the  jet  will  be  supersonic  and  this  problem  is  as  yet 
unsolved. 


The  approximate  differential  equation  used  is  (given  in  equa¬ 
tion  (84a)) 


8%  aq(l  +  bcr) 
(l  +  ca)^ 


(138) 


and  its  general  solution  has  been  found  to  be 


OO  ,  V 

i|r(cr,0)  =  (cxa  +  c2)(c3e  +  ^  (An  cos  n\0  +  Bn  sin  nA.©)wk^m  (z) 


n=l 


=  (cx0  +  c2)(c36  +  c^)  +  ^  Anwj^  (z)  sin  (nX.0  +  (%)  (139) 


n=l 


NACA  TN  2582 


66 


where  c^> 

c2’ 

c3> 

CV  V  BnJ 

,  and  X 

are 

constants;  An  Sin 

tf 

11 

> 

and  An  cos 

0^  = 

Bn 

are  obvious; 

2nX 
z  =  - 

c2 

<5b( 

1  +  cc)  ^  0, 

1 

1 

1 

P 

0 

1 

2b) 

li< 

0,  and  m  =  ^ 

1  (nX) 

^a(  c 

-  b)  , 

;  and  n  is  a 

posi- 

1 

U 

c^ 

tive  integer.  It  represents  the  required  solution  if  the  series 
converges. 


The  constants  cg,  c^,  c^,  Aq,  an,  and  X  are  determined 

from  the  boundary  conditions  which  are  shown  as  follows. 

60  arbitrary.-  The  flow  in  the  physical  plane  is  shown  in  fig¬ 
ure  16(a)  and  in  the  hodograph  planes  in  figures  16(b)  and  16(c).  Since 
it  is  known  that  the  discharge  Q  from  the  aperture  is  finite,  'the 
stream  function  t  is  bounded.  Moreover,  the  flow  is  symmetrical  with 
respect  to  the  center  line.  If  the  rate  of  total  discharge  is  introduced 
as  Q,  then  for  inside  the  vessel  there  can  be  written 


i  at  6  =  e0 

\|r  =  0  at  9  =  0  ► 

^  at  6  =  -0o 


(iko) 


Thus  it  is  obvious  that 
a  constant  q  and 


;  is  an  odd  function  of  0 .  Corresponding  to 
cr  on  the  outer  surface  of  the  jet  are  constant 


from  consideration  of  Bernoulli1 s  equation. 


Thus 


at 


a  =  cr0  =  Constant 
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First  of  all,  since  \|r  =  0  at  0  =  0, 


°n  =  0 


04  =  0 


Thus  equation  (139)  becomes 


(lte) 


00 

{cla  *  c2  )  cf  *  ]T 


Vm  <z> 


sin  nX6 


n=l 


(1^3) 


which  is  an  odd  function  of  0  as  required.  Next,  the  stream  function 
is  defined  in  the  interval  -0Q  %  0  <  0Q  or  the  period  20Q.  Thus 

(when  n  is  an  integer)  sin  nl(0  +  20o)  =  sin  nX0  or  2n\0o  =  2nit. 

*  =  r  d^) 

eo 

Substituting  equation  (l44)  into  equation  (1^3), 

00 

^(0,0)  =  (Cl«r  +  c2)c30  +  ^  AnW (z)  sin  |£ 0  (l45) 

n=i  0 


or 


Now,  since  \|r  is  bounded  everywhere,  particularly  when  a  — 
M  — ►  0,  must  be  equal  to  0.  Thus  equation  (l4-5)  yields 


t(a,0) 


=  Ao0  +  Anwk*i  sin  f- 

fcl  e° 


(1^6) 


0 
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or 


Q_  1 

mt  TT(n)  /  \ 

Wk,m  (zo) 


Substituting  into  equation  (lk8) } 


t(cr,0) 


W- 


(n) 

k,m 


nW- 


f(n) 

k,m 


(z) 

zo 


(151) 


Firsts  it  is  necessary  to  show  that  equation  (151)  is  valid  for  an 
incompressible  flow  which  can  be  considered  as  the  limiting  case  of 
a*  - ^00.  From  the  definition. 


a 


_P_  ia 
P0  a 


there  results,  for  incompressible  flow, 


a* 

lim  log£J  — 
a* — >00  4 


(152) 


and,  from  the  definition  of  z  for  very  large  values  of 


^since  =  lj 


=  2nkcr 


E 


From  equation  (87)  the  asymptotic  expansion  of  wj^  (z)  for  large 
values  of  a  is 


70 
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J* 


Introducing  the  above  results  into  equation  (153), 


Thus  the  stream  function  for  the  .incompressible  flow  is 


Q  JL 
2  e0 


(154) 


which  reduces  to  Chaplygin* s  case  if  0O 


is  set  equal  to  n/2. 


easily  proved,  because 
(z)  >  0  is  a  monotonically 

"k,m  \*o) 

decreasing  function  of  z  in  the  subsonic  range  (see  appendix  B). 

Compare  the  nth  term  of  equation  (151)  with  -  sin  nn  —  which  is  a 

n  9o 

term  of  a  convergent  series,  and  each  corresponding  term  is  smaller. 
Therefore,  equation  (l5l)  is  convergent. 


,(n) 

'k,m 


The  convergence  of  the  series  is 


Wil  1  (z) 


< 


(n) 


^  1  for 


\ 


z  \  zQ 


and  W 


(n) 

k,m 


The  main  interest  of  this  problem  has  been  the  minimum  width  of 
the  jet  which  occurs  at  X  =  00.  Therefore  the  solution  must  be  trans¬ 
formed  back  to  the  physical  plane. 

With  the  differential  relation  from  equations  (la)  and  (lb),  the 
following  equations  can  be  derived: 


M ,  Kft 

dcr  be 

hi 

be 


bo 


(155) 
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Substituting  the  above  equations  into  equation  ( 156)  , 


\K  A  1  ,  V~  Wk?n 
4  1  ®0  2  4r  M 


WW  <2>  ..  me 

/  /  \  cos  — — —  + 

-1 (-0)  s° 


2p  Q  v°°  w,  (n)(z) 

"i  0  \  nn0 

1  /  7 — 7 -  sm  -  dcr  + 

X  u(n)  /„  \  6, 


n=1  Wk,m  (Zo) 


J  2Q  V  k,: 
q  1  6o  ^  u(n 


r2-  «v(n)  U) 

V  — ,m  Bln  Hg  t 

n=l  wk,m  roj  0 


060  2  *  n=l  „<■>)  (z  ) 

L_  k,m  V  0/ 


W^.n  ^  ( z ) 

Jn)  /„  \  eo 


(158) 


which  is  the  conformal  relation  of  any  elementary  length  in  the  hodo- 
graph  plane  to  the  corresponding  element  in  the  physical  plane  as  long 

as  the  Jacobian  determinate  ji  0.  This  condition  is  automatically 

fuif ilied  in  the  subsonic  region  up  to  sonic  velocity  as  has  been  shown 
y  Tsien  (reference  45)  and  Craggs  (reference  48).  Since  dZ  is  an 
exact  differential,  there  can  "be  introduced, 


OK  1  ,  °l  <J> 

e°  2  “5=1  W<n)  (z.\ 


i2p  Q 


00  If 

yx 

-  T.r(n 


0  XL  (zo) 


(159a) 
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too  t  ( n)  /  x 

2QV~Wk,m  s.n  nnd  + 

do  n=l  (zQ)  6° 

k,m  V  0/ 


iPnQ1 


P0n  2 


s~~ 

1  .  > 

o  jL. —  __fnl 


Wv .m  nit 0 


=1  Wk,m  (zo) 


(159k) 


The  value  of  Z  can  he  determined  by  integrating  leaving  a  func 

tion  of  a  to  he  determined,  and,  consequently, 


Qe*y*f*w  (z)  f 
e°q  ^  «(!» h o ) 


(n)  t „\  HI  cos  -  i  sin  £6 


(€■ 


0OP1 


(n)  /_s  i  ^  sin  HI  6  -  cos  —  6 


Qpoe10  Wk,m  (  ^  eo  eo 


fci  WS  (zo) 


+  +  F(o) 


2©  0pq 


(160) 


where  F(<t)  is  to  he  determined.  This  can  he  done  hy  differentiating 
equation  (l60) 


SZ  =  Qje \ 

ha  Go*  fel 


fef  - (Zo) 


fXSrn  <z>  '  (z)] 


sin 


<z> 


)  (z)  +  £*  Kw[n>  (z)  l  +  K  +  F'  ( a) 

.  60  k>m  I  260q 


(161) 
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where 


Zrr  = 


2 rat 

' 


equations  (159a) 


and  equations  (l31a)  and  (l31b)  are  used, 
and  (161), 


Equating 


The  expressions  in  the  brackets  of  the  right-hand  side  of  equation  (l62) 
are  identically  equal  to  the  differential  equation;  equation  (78)  if 

o  n2  2 

V  =  — 5—  and  consequently  are  zero  as  shown. 

Thus;  choose 


F(ct)  =  Constant  =  Z 


(163) 


Now  equation  (l6o)  can  be  rewritten  in  a  more  suitable  form 

(g+i)e 


z  = 


-PpQ 

2p0oq 


16  +  JZ 

n=l 


0  p  *  (n)  ,  ,  TT(n)  il 

^  p0  Wk,m  (z)  +  ^k,m  (z)  e 


w[n)  (z  ) 
k,m  '  °/ 


—  +  1 


2  —  w3n^  (z)  -  (z) 

pQ  wk,m  wk,m  ^z'  g  \eQ  J 


"k°2  (zo) 


nn  -I 
eo  ' 


( l64a) 
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-PpQ 

2penq 


2£w* 


t!»  <*>  *  <*>  ^  (jr  +  i)e 


<-)  -  (*) 


(-0) 


(164c) 


TI16  convergence  of  "the  ahove  "two  series  must  he  established,  Firsts  ai 
cr  =  aQ  or  z  =  zQ}  equation  (l64b)  does  not  converge  for  6  =  0,  "because 

the  velocity  vector  of  q  =  q^  at  0=0  is  located  at  X  = 

At  very  large  values  of  n. 
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With  the  above  notations  the  asymptotic  expansion  of  -m  (z)  for 
large  values  of  n  can  be  introduced  as  follows: 


41 i  (z) 


h  r©r( 


k  + - m 

2 


a  - 1"0 


*  °(;E 


(166) 

at  sonic  velocity  (a  =  0),  For  the  subsonic  velocity  (cr  >  0),  a  similar 
expression  should  be  given;  the  case  a  =  0  is  more  critical.  Further¬ 
more  from  reference  6l,  page  352,  example  3> 


w'(n)  M 

k,m  v  1 
(z) 

k,m 


i)2l«inj  (»> 

J  k-l,m 


»t”i  <z> 


(167) 


Now  it  is  necessary  to  show  that,  for  large  values  of  n,  the  value  of 
w’(n)  (z) 

^,m -  does  not  depend  on  n  to  the  first  order.  Substituting 

W<”>  (z) 
k,m  v  7 

equation  (l66)  into  equation  (167), 


«k!m  (*)  g  1 

W<n)  (z)  2 

k,m 


E2  -  (*  - 1)2]  r(k  - 1  -  -) 


2  r(k  +  |  -  m) 

=  k  -  |  +  (m  +  k)e-45[l  +  0(|)] 


[ +  °(sl 


(168) 


which  is  a  constant  independent  of  l/n  to  the  first  order  for  large 
values  of  n.  Therefore  the  coefficients  involving  Whittaker* s  function 
for  large  values  of  n  can  be  written  as 


2p  Wk|m}  (zo)  +  ^2p[J 


P°  W(n)  ( z  \ 
Wk,m  \  o) 


r- 

p-i 


+  (m  +  k)e”^m  +  1 
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20  ^  (O 

P°WS(zo) 


-1~— k-i+(m  +  k)e“^“j 


(169b) 


both  of  which  are  bounded.  For  the  case  Qq  =  .  .  . ,  —  (s  is 

^  J  S 

an  integer)  the  series  in  equation  (l64c)  is  convergent.  This  can  be 
shown  at  the  surface  of  the  jet.  The  equation 


2«p  q 

OO  00 


sin  6  + 


p  u*  (  7  \ 

d  P0  Wk;m  \°J  n  sin  (sn  +  1  )G 
M  /  7  \  sn  +  1 


K  1  K) 

k,m  \  o  / 


a  Poo  „’(n)  , 

2  pQ  Wk,m  (zo) 

w[n)  (z  ) 
k,m  \  0/ 


sin  (sn  -  l)0 


sn  -  1 


(170) 


sp0Q  / 

is  dominated  by  -  1  +  E 

2itp  a  \ 


sin  n  0 


y2-  sin  nl©l\ 
n=l  n  / 


where  E 


and  F  are  positive  numbers  such  that  E  ^  2  -S2 


p.  W 


+  1  and 


0  C  M 


1 (n)  /  \ 

P  Wk.m  (zQ) 


F  >  0  —  k>m 

po  T.r(n) 


-  1  for  all  large  values  of  n, 


(zo) 


series 


Actually  fewer  terms  occur  in  equation  (170)  than  in  the  dominate 

V2"  Sil1  nl0l  1C  l0l  _ , _ ... _ /T7A\ 


V  sin  n 
/ —  n 


2  2 


Therefore  equations  (170)  and  conse¬ 


quently  equation  (l64c)  are  convergent  for  a  =  0.  Similarly,  the 
convergence  for  a  >  0  can  be  shown. 


11H 
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®o  ~  This  is  the  case  which  has  been  treated,  by  Chaplygin  and 

Lighthill.  The  condition  6  =  6Q  occurs  at  the  mouth  Of  the  aperture. 
Let  h/2  be  the  half  width  of  the  mouth.  Then, 


PpQ 


(-l)Dn 
kn2  -  1 


p.'m’  (  zo) 

Po'%i  (Zo) 


+ 


If,  at  X  =«,  the  width  of  the  jet  is  h^,  then,  by  definition, 

P  q 

^  =  (172) 

Mo 

Substituting  this  value  into  equation  (171), 


(173) 
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Making  use  of  the  result  of  the  well-known  Leibniz  series 
00  ^ 

)  ilil^ —  =  *  equation  (173)  can  he  written 
”1  2n  -  1  V 

t  2  (-Dn-  (z°>  *  £fe  _  A  +  £/i\ 

«  Pok^-l„(n)(2o)  ^  }  ”W 


l)nn  wk,n  (zo) 


,  1  +  16  p-  y  (-l)"n  »k,i 

*  pQ  Z_  4n2  -  1  u(n 


<nJ  (zo) 
k.m  \  0/ 


(nb) 


Just  as  a  check  for  the  case  of  incompressible  flow,  the  limiting  case 
of  a*  — ^  00  or  o'  — >  00  may  be  considered. 

First,  considering  equation  (l6T)  as  or  becomes  very  large  and  n 
is  finite. 


“k(,£)(*>  -  1 


-  1  "2  -  (k  -  I) 


2  z 


e  2 


e  2  zk 


k  -  i  -  m2  + 

2 


(e-^) 


(175) 


The  above  equation  tends  to  when  o  - — *  00,  owing  to  the  fact  that 

£  and  m  both  become  zero  as  a  — » 00  (z  — >00).  Substituting  the 
above  equation  into  equation  (17*0> 
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jl  =  1  +  a  f~  (-d11'1!! 
*  h°°  *  “ 1  ^n2  "  1 


=  1 


n=l 


2n’  -  1 


+  1 


(176a) 


Or,  inversely, 


00 

h~ 


2  +  it 


=  0.6110 


(176b) 


which  checks  with  the  well-known  result  of  Kirchoff. 

For  the  case  aQ  =  0  ( Mm  =  l)  the  width  ratio  h^/h  has  been 

calculated  from  equation  (17^).  It  is  found  to  be  0.7^6  by  summing 
terms  up  to  n  =  10.  It  is  expected  that  the  error  is  of  the 
order  -0.002.  Consequently,  it  checks  reasonably  well  with  the  Lighthill 
result  (hmjh  =  0.7^7)  which  is  based  on  the  exact  Chaplygin  function. 
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(«> 


Kummer  function 


p  1  n^X^a(c  -  b)  1  v/a  ab  .. 

m  =  - - -  ij  if  -  -  1 


—  _  m 
m  =  — 
z 


W<n)  U) 

k,m 


"><°)  <*> 
k,m  v  ' 

X,Y 

Z 


positive  integer 
pressure 

stagnation  pressure,  pressure  in  vessel 

free-stream  pressure,  pressure  surrounding  jet 

arbitrary  function  of  q  (Q(q)) 

rate  of  total  discharge 

local  velocity 

X  component  of  velocity 

Y  component  of  velocity 

Whittaker  function 


(n) 

'k,m 


dW^‘«  (z) 


dz 


coordinates  in  physical  plane 
complex  variable  (X  4-  iY) 


=  ~  ^ab(l  +  cff)  =  (1  +  c a)  if  =  1 


2n \ 
.  c 


ab 


2o  =  +  cao) 
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a,  |3 


constants 


3  = 


7  -  1 
7  =  Pp/py  =  4.4 


e  =  Q(q)  ^defined  by  e  =  - 


1  -  M‘ 


£  =  fi)  —  =  —  (supersonic) 

'2/  e3  q2 


-  /3\2  _  f£ 

~c»2 


-ers- 


X  =  -S- 


^0^0 


(subsonic) 


^1 


7+1 
7  -  1 


inclination  of  velocity  vector 

coordinate  of  source  location,  hodograph  plane 

slope  of  nozzle  wall 
characteristic  coordinate  (IT  +  0) 


coordinates  of  source  location  in  hodograph  plane 


characteristic  coordinate  (-IT  +  0) 


Po  1  -  M2 


dq 


positive  constant  dependent  on  boundary  conditions  (nX) 


local  density 
stagnation  density 
free- stream  density 


potential  function 
stream  function 


“  (£  -  i)l/g 


dq  (supersonic) 


Ua 


Q  + 


a/ 

V 


O  a* 


U* 


(l  -  m2)i/2  , 

-  dq 

q 


(subsonic) 


coordinate  of  source  location  in  hodograph  plane 


co  in  free  stream 


APPENDIX  B 


PROPERTIES  OF  WHITTAKER  FUNCTIONS 


Here  are  presented  a  few  important  properties  of  the  Whittaker 
function  which  are  useful  in  the  present  investigation.  Some  of  the 
properties  are  given  by  Sharma  in  reference  59*  He  shows  that: 

Theorem  I: 

The  functions  Wk^m  (z)  and  Wk^m+1  (z)  cannot  have 
a  common  zero  (root).  All  roots  of  Wk^m  (z)  are  simple; 

Wkm  (z)  and  Wk_-|^m  (z)  cannot  have  a  common  root. 

Between  any  two  consecutive  zeros  of  Wk^m  (z)  lies  one 
and  only  one  zero  of  Wk-1^m  (z). 

The  proof  is  rather  simple.  See  the  reference. 

Next,  consider  the  Whittaker  equation: 


(Bl) 


where 


W  =  W 


(n) 

k,m 


(z). 


the  Whittaker  function. 


In  the  present  case. 


(a  >  c  >  b  >  0)  =  1  .  (B2) 

c2 


where  n  is  a  positive  integer  and 


X 


is  a  positive  constant. 
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nr 


1  n2X.2/a  \  0,  p 

=  4  "  — (c  "  V  <  0  for  n  X  > 


4(a.  -  c) 


2  1  n2X.2/a 

nr-  =  -r-  -  -  1 

*  C2 


(^  -  1^  >  0  for  n2k2  < 


4(a  -  c) 


(B4) 


.2  1 


In  this  case  m  -  —  <  0;  the  powerful  Wintner  theorem  (refer- 
ence  74)  which  requires'  k  <  0  and.  m2  -  i  >  0  to  make 

4 

( 2 ) 

(“l)n  - -  £  0  (0  ^  z  ^  oo)  cannot  he  applied.  Actually, 


Wk^m  oscillate  in  the  range  -i  <  a  <  0.  But  the  monotonic 

(n) 

positive  nature  of  (z)  for  a  >  0  can  he  proved  as  follows: 

Theorem  II: 

For  nX  >  0  and  0^cx<oo  (l^M>0),  then 

wk^(z)>0- 

The  proof  is  as  follows: 

(\  2 

— j  =  4n2X2|— V  there  can  be 

■  .  \c  / 


written 


d^W  n^X^acr(l  +  ha) 
da2  (1  +  ca)2 


W  =  0 


or,  more  properly, 


.2,2 


dfW  =  n  X  acr(l  +  hor) 
da^  (l  +  ccr)2 


W 


(B5) 
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By  the  given  conditions,  as  long  as  cr  ^  0, 


i2\2aa(l  +  her)  > 


Therefore  — — S  and  W  have  the  same  sign.  But  W  —  (z)  is 

dex 


asymptotic  to  e' 2  zk  =  e'^0  En\(i  +  a)JnX(  c  +h)  >0  for  large 

values  of  a  and  tends  to  0  as  a  — »  0.  Hence  Wk^(z)  curves 

2_ 

downward  as  a  — >  oo  and  so  must  continue  to  do  so  while  W  — -  >  0. 

do 

This  proves  the  theorem. 

There  is  another  interesting  feature  of  the  Whittaker  equation. 

As  long  as  m2  is  real  and  the  function  representing  the  boundary  con- 
dition  is  real,  wj^  (z)  is  always  real,  no  matter  whether  m  is 
imaginary  or  not.  For  the  present  investigation,  m  is  imaginary,  hut 
wln^  (z)  is  real  (0  <  z  <  »).  This  question  has  puzzled  the  author 
for  some  time. 


>  0  for  large 
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Figure  2.-  Plot  of  K  as  a  function  of  t  for  0  <  t  <0.22  and 
log(-K)  as  a  function  of  t  for  0.22  <  t  <  1.0.  =  6  for 

7  =  1.4.  K  =  (%)2(1  -  M2)  =  — - M  =  0,.  K  =  1,  T  =  0 

(1  -  t)^1 

M  =  1,  K  =  0^  T  =  l/6  j  M  =  oo,  K  =  -oo^  T  =  1. 
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Figure  5.-  Plots  of  oo  and  Q  against  various  functions. 


NASA  TN  2582 


101 


.0123 


0  I  2  3  4  . 


M  — - 

(b)  As  functions  of  M.  See  equations  (l8c)  and  (22c).  M  =  0,  00  =  -ooj 

M  =  1,  cs  =  0;  M  =  1,  Cl  =  0;  M  =  Cl  =  2.277. 


Figure  5»-  Continued. 
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(c)  As  functions  of  a.  M  =  0,  a  -  ao,  00  =  M  =  1,  cr  =  0,  cu  =  0; 
M  =  1,  a  =  0,  ft  =  0;  M  =  00,  a  =  0;2513,  ft  =  2.277. 


Figure  5.-  Concluded. 


NACA  TN  2582 


Figure  8.-  Plot  of  ^  as  a  function  of  (jd  and  Q 
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Supersonic  approximations  in  comparison  with  exact 
1  against  ft*  Vertical  asymptote:  =  2.277* 
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cr  — 


0  .5  10  1.5  2.0 

^*■1-9)  (M».60)  (Mi. 29)  (M». 20) 


Figure  13.-  Detailed  comparison  of  third-order  approximation  with  exact 
K(a).  Solid  line, -exact;  triangle,  approximation.  Horizontal 
asymptote:  K  =  1;  vertical  asymptote:  a  =  -0.2513*  a  = 

b  =  1.61,  c  =  3.89.  2k  =  1. 

cd 


15H 


NACA  TN  2582 


111 


Figure  14.  -  Comparison  of  p0/p  in  hypothetical  gas  law  with  exact  gas 
law  as  a  function  of  ex.  Dotted  line,  exact ;  solid  line,  approximation; 
long  dashed  line,  Loewner  approximation  (reference  73).  Horizontal 

asymptote:  -£-  =  1.  —  .=  O.6339  at  a  =  0  (M  =  1). 
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Figure  15.-  Comparison  of  Mach  number  of  hypothetical  gas  and  exact  gas. 


Ob)  Hodograph  plane  (q  -  0).  (c)  Hodograph  plane  (a  -  0).  7  =  1 

*o  =  0;  =  0.5283;  ^  =  0.6339 

Po  Pc 

Figure  l6.~  Sketch  of  flow  through  an  inclined-walled  straight -edged 

aperture . 
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follow  the  exact  gas  law  closely  in  the  neighborhood  follow  the  exact  gas  law  closely  in  the  neighborhood 

of  sonic  velocity.  The  flow  through  an  aperture  is  of  sonic  velocity.  The  flow  through  an  aperture  is 

reexamined  by  comparing  the  present  approximate  reexamined  by  comparing  the  present  approximate 

relation  with  the  exact  Chaplygin  solutions-  relation  with  the  exact  Chaplygin  solutions- 
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